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Preface 

These notes are the second part of the tensor calculus documents which started with the 
previous set of introductory notes [11]. In the present text, we continue the discussion of 
selected topics of the subject at a higher level expanding, when necessary, some topics and 
developing further concepts and techniques. The purpose of the present text is to solidify, 
generalize, fill the gaps and make more rigorous what have been presented in the previous 
set of notes and to prepare the ground for the next set of notes. Unlike the previous 
notes which are largely based on a Cartesian approach, the present notes are essentially 
based on assuming an underlying general curvilinear coordinate system. We also provide a 
small sample of proofs to familiarize the reader with the tensor techniques inline with the 
tutorial nature of the present text; however, due to the limited objectives of the present 
text we do not provide comprehensive proofs and complete theoretical foundations for the 
provided materials. 

We generally follow the same conventions and notations used in the previous set of notes 
with the following amendments: 

• We use capital Gamma, T* fe , for the Christoffel symbols of the second kind which is 
more elegant and readable than the curly bracket notation |* fc } that we used in the 
previous notes insisting that, despite the suggestive appearance of the Gamma notation, 
the Christoffel symbols are not tensors in general. 

• Due to the restriction of using real (non-complex) quantities, as stated in the previous 
notes, all arguments of real-valued functions, like square roots and logarithmic functions, 
are assumed to be non-negative by taking the absolute value, if necessary, without using 
the absolute value symbol, as done by some authors. This is to simplify the notation and 
avoid confusion with the determinant notation. 

• We generalize the partial derivative notation so that <9* can symbolize the partial deriva- 
tive with respect to the u l coordinate of general curvilinear systems and not just for 
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Cartesian coordinates which are usually denoted by x l . The type of coordinates, being 
Cartesian or general or otherwise, will be determined by the context which should be 
obvious in all cases. 

• The summation symbol (i.e. Y2) is used in most cases when a summation is needed but 
the summation convention conditions do not apply or there is an ambiguity about it, e.g. 
when an index is repeated more than twice or a twice-repeated index is in an upper or 
lower state in both positions or a summation index is not repeated visually because it is 
part of a squared symbol. 

• “Tensor” and “Matrix” are not the same; however for ease of expression they are used 
sometimes interchangeably and hence some tensors may be referred to as matrices meaning 
the matrix representing the tensor. 

• In the present text, all coordinate transformations are assumed to be continuous, single 
valued and invertible. 
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1 Coordinate Systems, Spaces and Transformations 

• The focus of this section is coordinate systems, their types and transformations as well 
as some general properties of spaces which are needed for the development of the concepts 
and techniques of tensor calculus in the present and forthcoming notes. 

1.1 Coordinate Systems 

• In simple terms, a coordinate system is a mathematical device, essentially of geometric 
nature, used by an observer to identify the location of points and objects and describe 
events in generalized space which may include space-time. 

• The coordinates of a system can have the same or different physical dimensions. An 
example of the first is the Cartesian system where all the coordinates have the dimension 
of length, while examples of the second include the cylindrical and spherical systems where 
some coordinates have the dimension of length while others are dimensionless. 

• Generally, the physical dimensions of the components and basis vectors of the covariant 
and contravariant forms of a tensor are different. 

1.2 Spaces 

• A Riemannian space is a manifold characterized by the existing of a symmetric rank- 2 
tensor called the metric tensor. The components of this tensor, which can be in covariant 
(g t j) or contravariant ( g forms, are in general continuous variable functions of coordi- 
nates, i.e. g VJ = gijiu 1 , u 2 , . . . , u 11 ) and g l i = g l ^{u l , u 2 , . . . , u n ) where u 1 symbolize general 
coordinates. This tensor facilitates, among other things, the generalization of lengths and 
distances in general coordinates where the length of an element of arc, ds, is defined by: 

( ds ) 2 = g^dvldid 


( 1 ) 



1.2 Spaces 


In the special case of a Euclidean space coordinated by a rectangular system, the metric 
becomes the identity tensor, that is: 

9ii = 9 lJ = 9) = Sa = ^ = 5) (2) 

• The metric of a Riemannian space may be called the Riemannian metric. Similarly, the 
geometry of the space may be described as a Riemannian geometry. 

• All spaces dealt with in the present notes are Riemannian with well-defined metrics. 

• A manifold or space is dubbed “flat” when it is possible to find a coordinate system for the 
space with a diagonal metric tensor whose all diagonal elements are ±1; the space is called 
“curved” otherwise. Examples of flat space are the 3D Euclidean space coordinated by a 
rectangular Cartesian system whose metric tensor is diagonal with all the diagonal elements 
being +1, and the 4D Minkowski space-time whose metric is diagonal with elements of 
±1. Examples of curved space is the 4D space-time of general relativity in the presence 
of matter and energy. 

• When all the diagonal elements of the metric tensor of a flat space are +1, the space 
and the coordinate system may be described as homogeneous. 

• An nD manifold is Euclidean iff Rijki = 0 where Rijki is the Riemann tensor (see § 5.1); 
otherwise the manifold is curved to which the general Riemannian geometry applies. 

• A “held” is a function of the position vector over a region of space. Scalars, vectors 
and tensors may be dehned on a single point of the space or over an extended region 
of the space; in the latter case we have scalar fields, vector fields and tensor fields, e.g. 
temperature held, velocity held and stress held respectively. 

• In metric spaces, the physical quantities are independent of the form of description, being 
covariant or contravariant, as the metric tensor facilitates the transformation between the 
different forms; hence making the description objective. 
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1.3 Transformations 

• In general terms, a transformation from an nD space to another nD space is a corre- 
lation that maps a point from the first space (original) to a point in the second space 
(transformed) where each point in the original and transformed spaces is identified by n 
independent variables or coordinates. To distinguish between the two sets of coordinates 
in the two spaces, the coordinates of the points in the transformed space may be notated 
with barred symbols, e.g. (u 1 ,u 2 , . . . ,u n ) or («i,u 2 , • • • , h n ) where the superscripts and 
subscripts are indices, while the coordinates of the points in the original space are notated 
with unbarred similar symbols, e.g. (u 1 , u 2 , . . . , u n ) or (iq, n 2 , . . . ,u n ). Under certain 
conditions, which will be clarified later, such a transformation is unique and hence an 
inverse transformation from the transformed space to the original space is also defined. 
Mathematically, each one of the direct and inverse transformations can be regarded as 
a correlation expressed by a set of equations in which each coordinate in one space is 
considered as a function of the coordinates in the other space. Hence the transformations 
between the two sets of coordinates in the two spaces can by expressed mathematically by 
the following two sets of independent relations: 

u l = u\u l ,u 2 ,...,u n ) & = vi{u\u 2 ,...,u n ) (3) 

where i = 1,2, ... ,n with n being the space dimension. The independence of the above 
relations is guaranteed iff the Jacobian of the transformation does not vanish on any point 
in the space (see about Jacobian the forthcoming points). An alternative to viewing the 
transformation as a mapping between two different spaces is to view it as a correlation 
of the same point in the same space but observed from two different coordinate frames 
of reference which are subject to a similar transformation. The following points will be 
largely based on the latter view. 
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• As far as the notation is concerned, there is no fundamental difference between the 
barred and unbarred systems and hence the notation can be interchanged. 

• An injective transformation maps any two distinct points of the original space, ri and 
F 2 , onto two distinct points of the transformed space, fi and ? 2 - The image of an injective 
transformation, r, is regarded as coordinates for the point, and the collection of all such 
coordinates of the space points may be considered as a representation of a coordinate 
system for the space. If the mapping from an original rectangular system is linear, the 
coordinate system obtained from such a transformation is called “affine”. Coordinate 
systems which are not affine are described as “curvilinear” such as cylindrical and spherical 
systems. 

• The following n x n matrix of n 2 partial derivatives of the barred coordinates with 
respect to the unbarred coordinates is called the “Jacobian matrix” of the transformation 
between the barred and unbarred systems: 



du 1 

du 1 

du 1 

du 1 

du 2 

du n 

du 2 

du 2 

du 2 

du 1 

du 2 

du n 


du n 

du n 

du n 

du 1 

du 2 

du n 


( 4 ) 


while its determinant: 


J = det(J) 


( 5 ) 


is called the “Jacobian” of the transformation. 

• Barred and unbarred in the definition of Jacobian should be understood in a general 
sense not just as two labels since the Jacobian is not restricted to transformations between 
two systems of the same type but labeled as barred and unbarred. In fact the two coor- 
dinate systems can be fundamentally different in nature and not of the same type such 
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as Cartesian and general curvilinear. The Jacobian matrix and determinant represent 
any transformation by the above partial derivative matrix system between two coordi- 
nates defined by two different sets of coordinate variables not necessarily as barred and 
unbarred. The objective of defining the Jacobian as between barred and unbarred systems 
is generality and clarity. 

• The transformation from the unbarred coordinate system to the barred coordinate system 
is bijective 1 iff J 0 on any point in the transformed region of the space. In this case 
the inverse transformation from the barred to the unbarred system is also defined and 
bijective and is represented by the inverse of the Jacobian matrix: 2 3 

J = J - 1 ( 6 ) 


Consequently, the Jacobian of the inverse transformation, being the determinant of the 
inverse Jacobian matrix, is the reciprocal of the Jacobian of the original transformation: 

M (7) 

• A coordinate transformation is admissible iff the transformation is bijective with non- 
vanishing Jacobian and the transformation function is of class C 2 f 

• “Affine tensors” are tensors that correspond to admissible linear coordinate transforma- 
tions from an original rectangular system of coordinates. 

• Coordinate transformations are described as “proper” when they preserve the lianded- 


1 Bijective transformation means injective (one-to-one) and surjective (onto) mapping. 

2 Notationally, there is no fundamental difference between the barred and unbarred systems and hence 
the labeling is rather arbitrary and can be interchanged. Therefore, the Jacobian may be notated as 
barred over unbarred or the other way around. Yes, in a specific context when one of these is labeled 
as the Jacobian, the other one should be labeled as the inverse Jacobian to distinguish between the two 
opposite Jacobians and their corresponding transformations. 

3 C n continuity condition means that the function and all its first n partial derivatives do exist and 
are continuous in their domain of definition. Also, some authors impose a weaker condition of being of 
class C 1 . 
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ness (right- or left-handed) of the coordinate system and “improper” when they reverse 
the handedness. Improper transformations involve an odd number of coordinate axes 
inversions through the origin. 

• Inversion of axes may be called improper rotation while ordinary rotation is described 
as proper rotation. 

• Transformations of coordinates can be active, when they change the state of the observed 
object such as rotating the object in the space, or passive when they are based on keeping 
the state of the object and changing the state of the coordinate system which the object 
is observed from. In brief, the subject of an active transformation is the object while the 
subject of a passive transformation is the coordinate system. 

• An object that does not change by admissible coordinate transformations is described 
as “invariant” such as the value of a true scalar and the length of a vector. 

• As there are essentially two different types of basis vectors, namely tangent vectors of 
covariant nature and gradient vectors of contravariant nature, there are two main types 
of non-scalar tensors: contravariant and covariant tensors which are based on the type 
of the employed basis vectors. Tensors of mixed type employ in their definition mixed 
basis vectors of the opposite type to the corresponding indices of their components. As 
indicated earlier, the transformation between these different types is facilitated by the 
metric tensor. 

• A product or composition of coordinate transformations is a succession of transforma- 
tions where the output of one transformation is taken as the input to the next transfor- 
mation. In such cases, the Jacobian of the product is the product of the Jacobians of the 
individual transformations of which the product is made. 

• The collection of all admissible coordinate transformations with non-vanishing Jacobian 
form a group, that is they satisfy the properties of closure, associativity, identity and 
inverse. Hence, any convenient admissible coordinate system can be chosen as the point 
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of entry since other systems can be reached, if needed, through the set of admissible 
transformations. This is the cornerstone of building covariant physical theories which are 
independent of the subjective choice of coordinate systems and reference frames. 

• Transformation of coordinates is not a commutative operation. 

1.4 Coordinate Surfaces and Curves 

• The surfaces of constant coordinates at a certain point of the space meet to form curves 
(i.e. curves of intersection of these surfaces in pairs). The coordinate curves are these 
curves of mutual intersection of the surfaces of constant coordinates of the curvilinear 
system. 

• The above transformation equations (Eq. 3) are used to define the set of surfaces of 
constant coordinates and coordinate curves of mutual intersection of these surfaces. These 
coordinate surfaces and curves play a crucial role in the formulation and development of 
this subject. 

• The coordinate axes of a coordinate system can be rectilinear, and hence the coordinate 
curves are straight lines and the surfaces of constant coordinates are planes, as in the 
case of rectangular Cartesian systems, or curvilinear, and hence the coordinate curves are 
generalized curved paths and the surfaces of constant coordinates are generalized curved 
surfaces, as in the case of cylindrical and spherical systems. 

• In curvilinear coordinate systems, some or all of the coordinate surfaces are not planes 
and some or all of the coordinate lines are not straight lines. 

• Orthogonal coordinate systems are those for which the vectors tangent to the coordinate 
curves, as well as the vectors normal to the surfaces of constant coordinates, are mutually 
perpendicular at all points of the space. Consequently, in orthogonal coordinates, the coor- 
dinate surfaces are mutually perpendicular and the coordinate lines are also perpendicular 
at the point of intersection. 
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• In orthogonal coordinate systems, the corresponding covariant and contravariant basis 
vectors at any given point in the space are in the same direction, i.e. the tangent vector 
to a particular coordinate curve u l at a certain point and the gradient vector normal to 
the surface of constant u l at the same point have the same direction although they may 
be of different length. 

• A necessary and sufficient condition for a coordinate system to be orthogonal is that its 
metric tensor is diagonal. 

• An admissible coordinate transformation from a Cartesian system defines another Carte- 
sian system if the transformation is linear, and defines a curvilinear system if the trans- 
formation is nonlinear. 

1.5 Scale Factors 

• Scale factors (usually symbolized with hi, h 2 , ■ . . h n ) of a coordinate system are those 
factors which are required to multiply the coordinate differentials to obtain distances 
traversed during a change in the coordinate of that magnitude, e.g. p in the plane polar 
coordinate system which multiplies the differential of the polar angle dcj) to obtain the 
distance L traversed by a change of magnitude dcj) in the polar angle which is L — pdcj>. 
They are also used to normalize the basis vectors (refer to the previous notes [11] and 
forthcoming notes). 

• The scale factors for the Cartesian, cylindrical and spherical coordinate systems in 3D 
spaces are given in Table 1. 

• The scale factors are also used in the expressions for the differential elements of arc, 
surface and volume in general orthogonal coordinates, as described in § 2.3. 
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Table 1: The scale factors for the three most commonly used orthogonal coordinate systems 
in 3D spaces. The squares of these entries and the reciprocals of these squares give the 
diagonal elements of the covariant and contravariant metric tensors g tJ and g l] respectively 
of these systems. 



Cartesian (x, y, z ) 

Cylindrical (p, q b, z ) 

Spherical (r, 9, (ft) 

hi 

1 

1 

1 

h'2 

1 

P 

r 

h-.i 

1 

1 

r sind 


1.6 Basis Vectors in General Curvilinear Systems 

• The vectors providing the basis set for a coordinate system, which are not necessarily 
of unit length or mutually orthogonal, are of covariant type when they are tangent to 
the coordinate curves, and of contravariant type when they are perpendicular to the local 
surfaces of constant coordinates. Formally, the covariant and contravariant basis vectors 
are defined respectively by: 


Ei = 


dr 

dvi 


& E' : = V?/ 


( 8 ) 


where r is the position vector in Cartesian coordinates (a; 1 ,^ 2 , . . .), and u l are generalized 
curvilinear coordinates. 

• In general curvilinear coordinate systems, the covariant and contravariant basis sets, E, : 
and E', are functions of coordinates, i.e. 


E, : = Ej (u\ . . . , u n ) & E* = E' (w 1 , . . . , u n ) 


( 9 ) 


• Like other vectors, the covariant and contravariant basis vectors are related to each other 
through the metric tensor, that is: 


& E = g ij Ej 


Ej = <); , E 


( 10 ) 
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• The covariant and contravariant basis vectors are reciprocal basis systems, and hence 
in a 3D space with a right-handed coordinate system (ui,u 2 ,u 3 ) they are linked by the 
following relations: 


E 1 = 


Eo x E-., 


2 E 3 x Ei 


Ei • (E 2 x E 3 )’ 


E = 


3 El x E 2 


Ei • (Eo x E 3 


E = 


Ei • (Eo x E 3 


( 11 ) 


Ei = 


E 2 x E 3 


Eo = 


E 3 x E 1 


E, = 


E 1 x E 2 


E 1 • (E 2 x E 3 ) ’ E 1 • (E 2 x E 3 ) ’ ° E 1 • (E 2 x E 3 ) 

The relations in the last point may be expressed in a more compact form as follow: 


( 12 ) 


E' 


E; x Efc 
E • (Ej X Id;.. ) 


& 


E f 


E J x E fc 
E • (E i x E fc ) 


(13) 


where i,j, k take respectively the values 1, 2, 3 and the other two cyclic permutations (i.e. 
2, 3, 1 and 3, 1, 2). 

• The magnitude of the scalar triple product E* • (E ; x E&) represents the volume of the 
parallelepiped formed by E,, E ? and E/, : . 

• The magnitudes of the basis vectors in general orthogonal coordinates are given by: 


Ed — hi & 



(14) 


where hi is the scale factor for the i th coordinate. 

• The base vectors in the barred and unbarred general curvilinear coordinate systems are 
related by the following transformation rules: 


E * 

E' : 


dv? - 

& 


d 

dvi j 

Ej = 

dv? 3 

du l - • 

& 


du l ■ 

:E J 

dv? 

E' = 

E J 

du J 


( 15 ) 
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where the indexed u and u represent the coordinates in the nnbarred and barred sys- 
tems respectively. The transformation rules for the components can be straightforwardly 
concluded from the above rules; for example for a vector A which can be represented 
covariantly and contravariantly in the unbarred and barred systems as: 


A = E l Ai = E % 
A = E ; A = %A l 


(16) 


the transformation equations of its components between the two systems are given respec- 
tively by: 


A = 
A i = 


du J - 

& 


du 3 

du iAj 

Ai = 

du iAj 

du l - ■ 

& 


du 1 , ■ 

n 

du 3 

A 1 = 

n .A 3 
du 3 


(17) 


These transformation rules can be easily extended to higher rank tensors of different 
variance types, as detailed in the introductory notes [11]. 

• For a 3D manifold with a right-handed curvilinear coordinate system, we have: 


E, • (E 2 x E 3 ) = Jg & 


E 1 • (E 2 x E 3 ) = — 
1 ’ V9 


where g is the determinant of the covariant metric tensor, i.e. 


(18) 


g = fler !<),,) = 


(19) 


• Because E, • E y = g t] (Eq. 51) we have: 


J T J = [ 9 y] 


( 20 ) 


where J is the Jacobian matrix transforming between Cartesian and generalized coordi- 
nates, the superscript T represents matrix transposition, [g tJ \ is the matrix representing 
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the covariant metric tensor and the product on the left is a matrix product as defined in 
linear algebra which is equivalent to a dot product in tensor algebra. 

• Considering Eq. 20, the relation between the determinant of the metric tensor and the 
Jacobian is given by: 

g = J 2 (21) 

where J (= ||^ | with x for Cartesian and u for generalized coordinates) is the Jacobian 
of the transformation. 

• As explained earlier, in orthogonal coordinate systems the covariant and contravariant 
basis vectors, E, : and E', at any specific point of the space are in the same direction, and 
hence the normalization of each one of these basis sets, by dividing each basis vector by 
its magnitude, produces identical orthonormal basis sets . 4 * This, however, is not true in 
general curvilinear coordinates where each normalized basis set is different in general from 
the other. 

• When the covariant basis vectors E, are mutually orthogonal at all points of the space, 
we have: 

(A) the contravariant basis vectors E* are mutually orthogonal as well, 

(B) the covariant and contravariant metric tensors, g tJ and g l \ are diagonal with non- 
vanishing diagonal elements, i.e. 

9H = 9 ij = 0 j) ( 22 ) 

gn 7 ^ 0 & g 11 ^ 0 (no sum on i) (23) 

(C) the diagonal elements of the covariant and contravariant metric tensors are reciprocals, 

4 Consequently, there is no difference between the covariant and contravariant components of tensors 

with respect to such contravariant and covariant orthonormal basis sets. 
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g 11 = — (no summation) (24) 

Sii 

(D) the magnitude of the contravariant and covariant basis vectors are reciprocals, i.e. 

1 

r 

1.7 Covariant, Contravariant and Physical Representations 

• So far we are familiar with the covariant and contravariant (including mixed) repre- 
sentations of tensors. There is still another type of representation, that is the physical 
representation which is the common one in the applications of tensor calculus such as fluid 
and continuum mechanics. 

• The covariant and contravariant basis vectors, as well as the covariant and contravariant 
components of a vector, do not in general have the same physical dimensions as indicated 
earlier; moreover, the basis vectors may not have the same magnitude. This motivates 
the introduction of a more standard form of vectors by using physical components (which 
have the same dimensions) with normalized basis vectors (which are dimensionless with 
unit magnitude) where the metric tensor and the scale factors are employed to facilitate 
this process. The normalization of the basis vectors is done by dividing each vector by its 
magnitude. For example, the normalized covariant basis vectors of a general coordinate 
system, E*, are given by: 




E ; = 


E* 

IE; 


(no sum on i) 


(26) 


5 The comment “no summation” may not be needed in this type of expressions since both indices are 
of the same variance type in a generally non-Cartesian system. 
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which for an orthogonal coordinate system becomes: 


E, = 


E^_ 


Ej 

hi 


no snm on z 


(27) 


where gu is the i th diagonal element of the covariant metric tensor and hi is the scale factor 
of the i th coordinate as described previously. Consequently if the physical components of 
a vector are notated with a hat, then for an orthogonal system we have: 


A = A% = A% 



A 1 = ydfeA* = hiA 1 (no sum) (28) 


Similarly for the contravariant basis vectors we have: 


A i _ Ai 

\/fl i i hi 

(29) 

where g n is the i th diagonal element of the contravariant metric tensor. These definitions 
and processes can be easily extended to tensors of higher ranks. 

• The physical components of higher rank tensors are similarly defined as for rank-1 tensors 
by considering the basis vectors of the coordinated space where similar simplifications 
apply to orthogonal systems with mutually-perpendicular basis vectors. For example, for 
a rank- 2 tensor A with an orthogonal coordinate system, the physical components can be 
represented by: 

(no sum on i or j, with basis E'E 7 ) 

(no sum on i or j, with basis E,-Ey) (30) 

(no sum on i or j, with basis E,E J ) 


A-. - Aij 
hihj 


A 13 = hihjA 13 

Ai _ ^iA) 
A >~ — 


A = A;E' = A;E' : = Ai 




Ai = y^Ai 


On generalizing the above pattern, the physical components of a tensor of type (m, n) 
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in a general orthogonal coordinate system are given by: 

^ai...a m ^ a l ' ' ' ^ a m \ 

^■b! ...b n ~ h U ^b! ...b n 

n b\ ■ ■ • n b n 

• As a consequence of the last points, in a space with a well defined metric any tensor 
can be expressed in covariant or contravariant (including mixed) or physical forms using 
different sets of basis vectors. Moreover, these forms can be transformed from each other 
using the raising and lowering operators and scale factors. As before, for the Cartesian 
rectangular systems the covariant, contravariant and physical components are the same 
where the Kronecker delta is the metric tensor. 

• For orthogonal coordinate systems, the two sets of normalized covariant and contravari- 
ant basis vectors are identical as established earlier, and hence the physical components 
related to the covariant and contravariant components are identical as well. Consequently, 
for orthogonal systems with orthonormal basis vectors, the covariant, contravariant and 
physical components are identical. 

• The physical components of a tensor may be represented by the symbol of the tensor with 
subscripts denoting the coordinates of the employed coordinate system. For instance, if A 
is a vector in a 3D space with contravariant components A 1 or covariant components Ai, its 
physical components in Cartesian, cylindrical, spherical and general curvilinear systems 
may be denoted by (A x ,A y ,A z ), (A p ,A^,A z ), ( A r ,A e ,A and (A U ,A V ,A W ) respectively. 

• For consistency and dimensional homogeneity, the tensors in scientific applications are 
normally represented by their physical components with a set of normalized unit base vec- 
tors. The invariance of the tensor form then guarantees that the same tensor formulation 
is valid regardless of any particular coordinate system where standard tensor transforma- 
tions can be used to convert from one form to another without affecting the validity and 
invariance of the formulation. 
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2 Special Tensors 


• The subject of investigation of this section is those tensors that form an essential part of 
the tensor calculus theory, namely the Kronecker, the permutation and the metric tensors. 

2.1 Kronecker Tensor 

• This is a rank-2 symmetric, constant, isotropic tensor in all dimensions. 

• It is defined as: 



0 


1 


(* = 3) 


(* ^ j) 


(32) 


• The generalized Kronecker delta is defined as: 


( 


1 [(ji • • ■ j n ) is even permutation of [i\ . . . i n )] 

= — 1 [(j] . . .jn) is odd permutation of (R . . .i n )\ (33) 


0 


[repeated j’s] 


It can also be defined by the following n x n determinant: 



(34) 




where the <5* entries in the determinant are the normal Kronecker deltas as defined by Eq. 
32. 

• The relation between the rank-n permutation tensor and the generalized Kronecker delta 
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in an /rD space is given by: 




rl 2...n 


& 




°12...n 


(35) 


Hence, the permutation tensor e may be considered as a special case of the generalized Kro- 
necker delta. Consequently the permutation tensor can be written as an n x n determinant 
consisting of the normal Kronecker deltas. 

• If we define 


cij s;ijk 

° Im ~ ° Imk 


(36) 


then the well known e — 5 relation (Eq. 45) will take the following form: 



Win - FJi 


(37) 


Other identities involving 5 and e can also be formulated in terms of the generalized 
Kronecker delta. 


2.2 Permutation Tensor 

• This tensor has a rank equal to the number of dimensions of the space. Hence, a rank-n 
permutation tensor has n n components. 

• ft is a relative tensor of weight —1 for its covariant form and +1 for its contravariant 
form. 

• ft is isotropic and totally anti-symmetric in each pair of its indices, i.e. it changes sign 
on swapping any two of its indices. 

• ft is a pseudo tensor since it acquires a minus sign under improper orthogonal transfor- 


mation of coordinates. 
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The rank-n permutation tensor is defined as: 


lll2 ■■■In / 

t — c UZ 2 ...in 


1 [(ii, i 2 , . . . , i n ) is even permutation of (1,2,..., n)] 

— 1 [(ii, i 2 , ■ ■ ■ , i n ) is odd permutation of (1,2,..., n)] 

0 [repeated index] 


(38) 


For the rank-n permutation tensor we have: 


.6 


v 'CL\a2--CLn 


n— 1 

n 

i=l L j=i~\~ 1 


Z! 


J J i a j a i ) 


5(n - 1) 


[| (%■ - a i ) 


(39) 




where 5(n — 1) is the super-factorial function of (n — 1) which is defined by: 


k 

S{k) = IJi! = l!-2!-...-fc! 

2=1 


(40) 


• A simpler formula for the rank-n permutation tensor can be obtained from the previous 
one by ignoring the magnitude of the multiplication factors and taking their signs only, 
that is: 


where 


£a\<i2---a n 

= n - 

l<2<j<n 

y CLj Cli) 

(41) 


( 

+i 

(. k > 0) 


cr(k) = < 

-l 

(k < 0) 

(42) 


0 

(* = o) 





• The sign function in the previous point can be expressed in a more direct form by dividing 
each argument of the multiplicative factors in Eq. 41 by its absolute value, noting that 


6 The following formulae also apply to the contravariant form. 
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none of these factors is zero, and hence Eq. 41 becomes: 


^aiCL2---a n 


n 


( CLj &i) 


\CLj — CLi 
l<i<j<n ^ 


(43) 


For the rank-3 permutation tensor we have: 


e ijk e lmn = 


5* 

°l °m °n 


si si, si 


r) k r) k 
°l °m °n 


(44) 


c ijk 


= S‘A - S‘ m Sj 


(45) 


For the rank-n permutation tensor we have: 


e, lh ..., n = n\ 


(46) 


For the rank-n permutation tensor we have: 




3132 "On 




• S ) 1 

Jn 

• 

Jn 


(47) 
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6? 

Jn 


• On comparing Eqs. 34 and 47 we obtain the following identity: 


...t n Jll •••In ^ 

U jl-jn ~ L ^jl-jn 


(48) 


• Based on the previous point, the generalized Kronecker delta is the result of multiplying 
two relative tensors one of weight w — +1 and the other of weight w = — 1 and hence the 
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generalized Kronecker delta has a weight of w = 0; therefore the generalized Kronecker 
delta is an absolute tensor. ' 

• As has been stated previously, e ll '" ln and are relative tensors of weight +1 and 

— 1 respectively. It is desirable to define absolute covariant and contravariant forms of the 
permutation tensor, marked with underline, by the following relations: 8 

<L h ...i n vVe.-.t. & = _L e u...i„ (49) 

where g is the determinant of the covariant metric tensor g pq . 

• The e — S identity (Eqs. 37 and 45) can be generalized by employing the metric tensor 
with the absolute permutation tensor: 

9 ^ 9ikl9jmn — 9km9ln 9kn9lm (50) 


2.3 Metric Tensor 

• One of the main objectives of the metric, which is a rank-2 symmetric absolute non- 
singular 9 tensor, is to generalize the concept of distance to general curvilinear coordinate 
frames and hence maintain the invariance of distance in different coordinate systems. This 
tensor is also used to raise and lower indices and thus facilitate the transformation between 
the covariant and contravariant types. 

• In general, the coordinate system and the space metric are independent entities. Yes, 
some coordinate systems may be defined by having a specific metric in which case the two 


'The multiplication of relative tensors produces a tensor whose weight is the sum of the weights of 
the original tensors. 

8 The contravariant form requires a sign function with details out of scope of the present text (see 
[13]); however, for the rank-3 permutation tensor which is the one used mostly in the forthcoming notes 
the above expression stands as it is. 

9 We mean that the matrix representing the tensor is invertible and hence its determinant does not 
vanish at any point of the space. 
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are correlated. This is the case in the Cartesian coordinate systems which are based in 
their definition on presuming an underlying Euclidean metric. 

• The components of the metric tensor are given by: 

g l3 = Ej • E j & g ij = E* • E J (51) 


where the indexed E are the covariant and contravariant basis vectors as defined previously 
in § 1.6. Because of these relations, the vectors Ej and E' may be denoted by g,; and g' 
respectively which is more suggestive of their relation to the metric tensor. 

• As a consequence of the last point, the covariant metric tensor can also be defined as: 


dx k dx k 
du 1 dui 


(52) 


where 

x k = x k (u 1 , . . . , u 11 ) (k — 1, ...,n) (53) 

are independent coordinates in an nD space with a rectangular Cartesian system, and 
u l (i — l,...,n) are independent generalized curvilinear coordinates. Similarly for the 
contravariant metric tensor we have: 


d u l du j 
dx k dx k 


(54) 


• The coefficients of the metric tensor may also be considered as the components of the 
unit tensor in its two variance forms, that is: 


5 = ()jj E'E' = //TIE 


(55) 


As stated already, the basis vectors, whether covariant or contravariant, in general 
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coordinate systems are not necessarily mutually orthogonal and hence the metric tensor 
is not diagonal in general since the dot products given in Eqs. 51, 52 and 54 are not 
necessarily zero when i ^ j. Moreover, since those basis vectors are not necessarily of unit 
length, the entries of the metric tensor are not of unit magnitude in general. However, 
since the dot product of vectors is a commutative operation, the metric tensor is necessarily 
symmetric. 

• The entries of the metric tensor, including the diagonal elements, can be positive or 
negative. 

• The covariant and contravariant forms of the metric tensor are inverses of each other 
and hence: 

g ik g kj = S‘ j k g ik g k > = S’ (56) 

where these equations can be seen as a matrix multiplication (row x column). 

• A result from the previous points is that: 


(E* • E') (E, • Efc) = g*g jk = S\ 

(57) 

(E* • Ej) (E' • E k ) = g l3 (f k = 8 k 

• As the metric tensor has an inverse, it should not be singular and hence its determinant, 
which is in general a function of coordinates like the metric tensor itself, should not vanish 
at any point in the space, that is: 


giu 1 ,.. .,«”) = det (g y ) / 0 


(58) 


• The mixed type metric tensor is given by: 


q l 3 = E* • Ej = & 


3 — 


= E; • W = 8- 


(59) 
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and hence it is the identity tensor. These equations represent the fact that the covariant 
and contravariant basis vectors are reciprocal sets. 

• From the previous points, it can be concluded that the metric tensor is in fact a trans- 
formation of the Kronecker delta in its different variance types from a rectangular system 
to a general curvilinear system, that is: 


dx k dx l dx k dx k 

du * dip kl du l dip 1 j 

q ij diP_diP_ E?; EJ 

dx k dx l dx k dx k 

i = dv^dx[ §k = dv^dx^ = e* • E- 
3 dx k dui 1 dx k dui J 


( co variant) 

(contravariant) 

(mixed) 


• Because of the relations: 


( 60 ) 


A‘ = A ■ E‘ = AjW ■ E' = 

At = A ■ E f = A’Ej ■ E ; = tig# 


( 61 ) 


the metric tensor is used as an operator for raising and lowering indices and hence facili- 
tating the transformation between the covariant and contravariant types of vectors. By a 
similar argument, the above can be easily generalized where the contravariant metric ten- 
sor is used for raising covariant indices and the covariant metric tensor is used for lowering 
contravariant indices of tensors of any rank, e.g. 


A\ = g^Ajk k A, a = 


( 62 ) 


Consequently, any tensor in a Riemannian space with well-defined metric can be cast into 
covariant or contravariant or mixed forms . 10 

• In the raising and lowering of index operations the metric tensor acts, like a Kronecker 
delta, as an index replacement operator as well as shifting the index position. 


10 For mixed form the rank should be > 1. 
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• In general, the order of the raised and lowered indices is important and hence 

g ik A jk = A/ and g ik A kj = A i j (63) 

are different unless the tensor is symmetric in its two indices, i.e. Aj k = A k j. A dot may 
be used to indicate the original position of the shifted index and hence the order of the 
indices is recorded, e.g. A,- 1 , and A.® • for the above examples respectively, although this is 
redundant in the case of symmetry. 11 

• Raising and lowering of indices is a reversible process; hence keeping a record of the 
original position of the shifted indices will facilitate the reversal. 

• For a space with a coordinate system in which the metric tensor can be cast into a 
diagonal form with all the diagonal entries being of unity magnitude (i.e. ±1) the metric 
is called flat. 

• If g and g are the determinants of the covariant metric tensor in the unbarred and barred 
systems respectively, i.e. g = det ( gij ) and g = det (fjij), then 

g = J 2 g & VH=JVg (64) 

where J (= |§||) is the Jacobian of the transformation between the unbarred and barred 
systems. Consequently, the determinant of the covariant metric and its square root are 
relative scalar invariants of weight +2 and +1 respectively. 

• A “conjugate” or “associated” tensor of a tensor in a metric space is a tensor obtained 
by inner product multiplication, once or more, of the original tensor by the covariant or 
contravariant forms of the metric tensor. 

• All tensors associated with a particular tensor through the metric tensor represent the 

11 Dots may also be inserted in the tensor symbols to remove any ambiguity about the order of the 
indices even without the action of the raising and lowering operators. 
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same tensor but in different reference frames since the association is no more than raising 
or lowering indices by the metric tensor which is equivalent to a representation of the 
components of the tensor relative to different basis sets. 

• A sufficient and necessary condition for the components of the metric tensor to be 
constants in a given coordinate system is that the Christoffel symbols of the first or second 
kind vanish identically (refer to 3.1). 

• The metric tensor behaves as a constant with respect to covariant and absolute differen- 
tiation (see § 3.2 and § 3.3). Hence, in all coordinate systems the covariant and absolute 
derivatives of the metric tensor are zero; moreover, the covariant and absolute derivative 
operators bypass the metric tensor in differentiating inner and outer products of tensors 
involving the metric tensor. 

• In general orthogonal coordinate systems in nD spaces the metric tensor and its inverse 
are diagonal, that is: 


9a = g'’ = o (i f- j) 


(65) 


moreover, we have: 



no sum on i 


( 66 ) 


det (g is ) = g = 9n922 ■ ■ 9nn - ff (hif 


(67) 


-l -1 


det (g lJ ) 


1 


1 


II (^) 2 


( 68 ) 


9 911922 ■ ■ ■ 9 Tin 


where hi (= |E ( -|) are the scale factors, as described previously. 


• A Riemannian metric, r/., ? . in a particular coordinate system is a Euclidean metric if it 
can be transformed to the identity tensor, <5^, by a permissible coordinate transformation. 

• The Minkowski metric, which is the metric tensor of special relativity, is given by one 
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of the following two forms: 


[ j «] = [g'’_ 


1 

0 

0 

1 — 

o 


i 

i 

0 

0 

1 

o 

0 

1 

0 

0 

\9a\ = [g 13 ] = 

0 

1 

0 

0 

0 

0 

-1 

0 


0 

0 

1 

0 

0 

0 

0 

-1 


0 

0 

0 

1 


(69) 


Consequently, the line element ds can be imaginary. 

• The partial derivatives of the covariant and contravariant metric tensors satisfy the 
following identities: 


dk9ij flm.j finish 9 


(70) 


dj 1 = -g mi g ,n 8 h g nm 

• In the following subsections, we investigate a number of mathematical objects whose 
definitions and applications are dependent on the metric tensor. 


2.3.1 Dot Product 

• The dot product of two basis vectors in general curvilinear coordinates was given earlier 
in this section. This will be used in the following points to develop expressions for the dot 
product of vectors and tensors in general. 

• The dot product of two vectors, A and B, in general curvilinear coordinates using their 
covariant and contravariant forms, as well as opposite forms, is given by: 


A B = AiK • BjEP = AiBj E i ■ E j = g lJ A l B J = A j P> l = A, PA 

A B = A% ■ B% = .T/CK, • Ej = g ij A i B j = AjB j = A' B>, (71) 

A B = A Al ■ B J Ej = A i B :i E i • E, = S)AiB j = AjB j 

= .TIC • BjEP = A i BjE i ■ E = h/.l'/i, = A l B t 


A B 
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In brief, the dot product of two vectors is the dot product of their two basis vectors 
multiplied algebraically by the algebraic product of their components. Because the dot 
product of basis vectors is a metric tensor, the metric tensor will act on the components by 
raising or lowering the index of one component or by replacing the index of a component. 
• The dot product operations outlined in the previous point can be easily extended to 
tensors of higher ranks where the covariant and contravariant forms of the components 
and basis vectors are treated in a similar manner to the above to obtain the dot product. 
For instance, the dot product of a rank-2 tensor of contravariant components A ^ and a 
vector of covariant components Bk is given by: 


A B = (A ij E,Ej) ■ (B k E k ) = A ij B k (Ivlv • E fc ) = A ij B k E,S- = A ; 'P> ; K: (72) 

that is, the i th component of this product, which is a contravariant vector, is: 

[A • Bf = A ij Bj (73) 

• From the previous points, the dot product in general curvilinear coordinates occurs 
between two vectors of opposite variance type. Therefore, to obtain the dot product of 
two vectors of the same variance type, one of the vectors should be converted to the 
opposite type by the raising/lowering operator, followed by the inner product operation. 
This can be generalized to the dot product of higher-rank tensors where the two contracted 
indices of the dot product should be of opposite variance type and hence the index-shifting 
operator in the form of the metric tensor should be used, if necessary, to achieve this. 

• The generalized dot product of two tensors is an invariant under permissible coordinate 


transformations. 
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2.3.2 Cross Product 

• The cross product of two covariant basis vectors in general curvilinear coordinates is 
given by: 


^ ^ dx l dx m 

E;XE ^=^ e ' X 9^ 


Gm — 


dx l dx m 
du i dui 


6/ X 0 m 


dx l Q x ™ 

du 1 dui 


^Imn^n 


(74) 


where the indexed x and u are the coordinates of Cartesian and general curvilinear sys- 
tems respectively, the indexed e are the Cartesian base vectors 12 and ei mn = e lmn is the 
permutation relative tensor as defined in Eq. 38. Now since e n = e n = %^E k , the last 


equation becomes: 


^ ^ dx l dx m dx 11 _ k _ k 

E,- x E, = ——-—ei mn E k = e iik E k 


^ du* dui du k " Lmn ^ - ijk 
where the underlined absolute covariant permutation tensor is defined as: 


(75) 


dx 1 Q x m Q x n 

-‘ ik = ThX 


(76) 


So the final result is: 

E. x E, = e m E k (77) 

By a similar reasoning, we obtain the following expression for the cross product of two 
contravariant basis vectors in general curvilinear coordinates: 


E' x E' = e«‘Ej 


(78) 


where the absolute contravariant permutation tensor is defined by: 


e ijk = ^ duj ^ c lmn 

dx 1 dx m dx n 

12 For Cartesian systems, there is no difference between covariant and contravariant tensors and hence 
e, : = e*. We also note that for Cartesian systems g = 1. 
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• Considering Eq. 49, the above equations can also be expressed as: 

E i x E j = e ijk E k = y/_ ge ijk E k 

ijk 

E' x E 3 = e l3k E k = — E k 

V9 

where e l]k = e l3k are as defined previously (Eq. 38). 

• The cross product of non-basis vectors follows similar rules to those outlined above for 
the basis vectors; the only difference is that the algebraic product of the components is 
used as a scale factor for the cross product of their basis vectors. For example, the cross 
product of two contravariant vectors, A 1 and B 3 , is given by: 

A x B = (A*Ej) x (B'Ej) = (E ; x E ? ) = e ijk A i B- i E k (82) 

that is, the k th component of this product, which is a vector with covariant components, 
is: 

[A x B]* = e ijk A‘B‘ (83) 

Similarly, the cross product of two covariant vectors, Ai and Bj, is given by: 

A x B = (A,E' : ) x (B :j E‘) = A,Bj (E f x E J ) = e lik A i B J E t (84) 

with the k th contravariant component being given by: 


(80) 


(81) 


[A x B]* = e vt A t Bj 


(85) 
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2.3.3 Line Element 

• The displacement differential vector in general curvilinear coordinate systems is given 
by: 

dr = ~z—du l = E idu 1 = E, : | — 9-du l = E, : E idu 1 ( 86 ) 

C / LL J—i 2 

l 11 l 

where r is the position vector as defined previously. 

• The line element ds, which may also be called the differential of arc length, in general 
curvilinear coordinate systems is given by: 

(ds) 2 = dr ■ dr = E idu 1 ■ E jdv? = (E, • E y ) du l du = g^diddu 3 (87) 


where g l3 is the covariant metric tensor. 

• For orthogonal coordinate systems, the metric tensor is given by: 


9ij 


0 



(* + j) 
(* = j) 


( 88 ) 


where hi is the scale factor of the respective coordinate u l . Hence, the last part of Eq. 87 
becomes: 

(ds) 2 = (hi) 2 du l du l (89) 

i 

with no cross terms (i.e. terms of products involving more than one coordinate like du l du i 
where i ^ j) which are generally present in the case of non-orthogonal curvilinear systems. 
• On conducting a transformation from one coordinate system to another coordinate 
system, marked with barred coordinates, u, the line element will be expressed in the new 
system as: 

(ds) 2 = cjijdu l du 


(90) 
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Since the line element is an invariant quantity, the same symbol ( ds ) 2 is used in both Eqs. 
87 and 90. 


2.3.4 Surface Element 

• In general curvilinear coordinates of a 3D space, an infinitesimal element of area on the 
surface u 1 = ci, where ci is a constant, is obtained by taking the magnitude of the cross 
product of the displacement vectors in the directions of the other two coordinates on that 
surface. Hence, the generalized differential of area element on the surface u 1 = C\ is given 
by: 

dA{u l = Ci) = \dr 2 x dr 3 | 

= IE-2 x E 3 | du 2 du 3 


2 •*> I-/3 | 

= |e 231 E 1 | du 2 du 3 


= |e 231 | E 1 1 du 2 du 3 


(Eq. 77) 


(91) 


(Eqs. 49 & 98) 
(Eq. 51) 


= ^/gV E 1 • E 1 du 2 du 3 

— yjg du 2 du 3 

— yj gg 1 1 du 2 du 3 

• On generalizing the above argument, the differential area element in a 3D space on the 
surface u l — c t (i — 1, 2, 3) where c* is a constant is given by: 


dA{u l = Ci) = \J gg n du J d;u k 


{i 7 ^ j 7 ^ k, no sum on 1 ) 


(92) 


In general orthogonal coordinates in a 3D space we have: 


= J (h i ) 2 (/i J ) 2 (h fe ) 2 7 — = hjh k 


(i / j ^ k, no sum on any index) (93) 
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and hence Eq. 92 becomes: 


dA(u l = Cj) = hjhkdu^ du k 


k, no sum on any index) (94) 


The last formula represents the area of a surface differential with sides h 3 d:u? and hkdu k 
(no sum on j, k ). 


2.3.5 Volume Element 

• In general curvilinear coordinates of a 3D space, an infinitesimal element of volume, 
represented by a parallelepiped spanned by the three displacement vectors dr^ (i = 1, 2, 3), 
is obtained by taking the magnitude of the scalar triple product of these vectors. Hence, 
the generalized differential volume element is given by: 


|dri ■ (dr 2 x dr 3 )| 


Ex • (E 2 x E 3 ) du l du 2 du 3 


Ei • e^iE 1 dv}du 2 du 3 

(Eq. 77) 

Ei • E 1 e 23 i du l du 2 du 3 


di 1 e 23 i \dv}du 2 dv? 

(Eq. 59) 

sj g du l du 2 du z 

(Eq. 49) 

J dv}du 2 du 3 

(Eq. 21) 


(95) 


where g is the determinant of the covariant metric tensor < 7 ^, and J is the Jacobian 13 of 
the transformation as defined previously. The last line in the last equation is particularly 
relevant to the case of change of variables in multivariate integrals where the Jacobian 
facilitates the transformation. 

1 ! Due to the freedom of choice in the order of the variables, which is related to the choice of the system 

handedness hence affecting the sign of the determinant Jacobian, the sign of the determinant should be 

adjusted if necessary to have a proper sign for the volume element. 
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• The formulae in the last point for a 3D space can be extended to the differential of 
generalized volume element 14 in general curvilinear coordinates in an nD space as follow: 

dV = y/gdu 1 . . . du n = J du 1 . . . du n (96) 

• In general orthogonal coordinate systems in a 3D space, the above formulae become: 

dV = h\h 2 h 3 du 1 du 2 du 3 (97) 

where hi,h 2 and h 3 are the scale factors. The last formula represents the volume of a 
parallelepiped with edges hidu 1 , h 2 du 2 and h 3 du 3 . 

2.3.6 Magnitude of Vector 

• The magnitude of a contravariant vector A is given by: 

|A| = y/A ■ A = ^(E, • Ej) A l Ai = sJg^Ao = yj AjAi = y/A*Ai (98) 

A similar expression can be obtained for the covariant form of the vector, that is: 

| A | = V A ■ A = J (E' • E-i) AiAj = yJlfUgAj = \/A j Aj = \[AJd (99) 

The magnitude of a vector can also be obtained more directly from the dot product of the 
covariant and contravariant forms of the vector: 

|A| = Va • A = yj (E' : • Ej) AiAi = sj^A^ = y/ A i A i = ^AjAi (100) 


14 Generalized volume elements are used, for instance, to represent the change of variables in multi- 
variable integrations. 
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2.3.7 Angle Between Vectors 

• The angle 9 between two contravariant or two covariant vectors A and B is given 
respectively by: 


A B gijA'B 3 g lJ AiBj 

cos u — — = ^ _ 

|A||B| a / g k iA k A l yj g mn B m B n \J g kl A k A ly J g mn B m B n 

For two vectors of opposite variance type we have: 


( 101 ) 


A B _ r _ \jr 

A||B| ^/,j M A k A'^(f"'B m B n y g kl A k A lx /g mn B m B’‘ 


( 102 ) 


2.3.8 Length of Curve 

• In general curvilinear coordinates, the length of a f-parameterized space curve r(£) 
dehned by u 1 = u l (t), which represents the distance traversed along the curve on moving 
between its start point S and end point E , is given by: 15 


L = 



a J gijdvtdvj 



du l dui 

gii HH dt 


(103) 


where £ is a scalar variable parameter, and £i and t 2 are the values of £ corresponding to 
the start and end points respectively. 

• The length of curve is used to define the geodesic which is the path of the shortest distance 
connecting two points in a Riemannian space. Although the geodesic is a straight line in 
a Euclidean space, it is a generalized curved path in a general Riemannian space. 


15 Some authors add a sign indicator to ensure that the argument of the square root is positive. However, 
as indicated in the Preface, such a condition is assumed when needed since we deal with non-complex 
values only. 
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3 Covariant and Absolute Differentiation 

• The focus of this section is the investigation of covariant and absolute differentiation 
operations which are closely linked. These operations represent generalization of tensor 
differentiation in general curvilinear coordinate systems. Briefly, the differential change 
of a tensor in general curvilinear coordinate systems is the result of a change in the 
base vectors and a change in the tensor components. Hence, covariant and absolute 
differentiation, in place of the normal differentiation, are defined and employed to account 
for both of these changes. Since Christoffel symbols are crucial in the formulation and 
application of covariant and absolute differentiation, the first subsection of the present 
section will be dedicated to these symbols and their properties. 

3.1 Christoffel Symbols 

• We start by investigating the main properties of the Christoffel symbols which play 
crucial roles in tensor calculus in general and are needed for the subsequent development 
of the present and forthcoming sections as well as the future notes. 

• Christoffel symbols are classified as those of the first kind and those of the second kind. 
These two kinds are linked through the index raising and lowering operators. Both kinds 
of Christoffel symbols are variable functions of coordinates in general. 

• Christoffel symbols of the first and second kind are not tensors in general although they 
are affine tensors of rank-3. 

• As a consequence of the last point, if all the Christoffel symbols of either kind vanished 
in a particular coordinate system they will not necessarily vanish in other systems; for 
instance they all vanish in Cartesian systems but not in cylindrical or spherical systems, 
as has been established previously [11] and will be investigated further in the forthcoming 
points. 
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• Christoffel symbols of the first kind are given by: 

[ij, l \ = \ ( d i9u + diQji ~ diQij ) (104) 

where the indexed g is the covariant form of the metric tensor. 

• Christoffel symbols of the second kind are obtained by raising the third index of the 
Christoffel symbols of the first kind, that is: 

kl 

r % = g kl [? ij , l \ = Y (®j9u + d%gn - digij ) (105) 

where the indexed g is the metric tensor in its contravariant and covariant forms with 
implied summation over /. 

• Similarly, the Christoffel symbols of the first kind can be obtained from the Christoffel 
symbols of the second kind by reversing the above process through lowering the upper 
index, that is: 

Skmltf = g km g kl [ij, l ] = 8 l m [ij, l ] = [ij, m] (106) 

• For an nD space with n covariant basis vectors (Ei, E 2 , . . . , E„) spanning the space, the 
derivative <9jE, for any given i is a vector within the space and hence it is in general a 
linear combination of all the basis vectors. The Christoffel symbols of the second kind are 
the components of this linear combination, that is: 

djEi = rf.E fc (107) 

Similarly, for the contravariant basis vectors we have: 

djE' = - 1 '[■ IT 


(108) 
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• By inner product multiplication of the previous relations with the basis vectors we 
obtain: 

E' ■ d k E j = V) k & E, ■ d k E j = -rj k (109) 

Similarly: 

Efc • djE, = g mk E m • djEi = g mk T ™ = [ij, k] (110) 

• Christoffel symbols of the first and second kind are symmetric in their paired indices, 
that is: 

[%j,k] = \n,k] k r*=r‘ (in) 

• The partial derivative of the components of the covariant metric tensor and the Christof- 
fel symbols of the first kind satisfy the following identity, which is essentially based on the 
forthcoming Ricci Theorem: 

djgu = [ij,l} + [jl,i\ (112) 

This relation can also be written in terms of the Christoffel symbols of the second kind 
using the index shifting operator: 


d]fhi — 9ki^ij + 9ki^ k ji (113) 

A related formula for the partial derivative of the components of the contravariant metric 
tensor, which can be obtained by partial differentiation of the relation 9i m g m J = Sj with 
respect to the k th coordinate, is given by: 

9imd k 9 m3 = - 9 mj dkgim (114) 


• Christoffel symbols of the second kind with two identical indices of opposite variance 
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type satisfy the following relations: 


= I# = Y g dig = 2 9i ^ hl 9 ^ = &i ^ hl ^ = ~^ di ^ 


1 

V9' 


(115) 


where the main relation can be derived as follow: 


Q jl 

T ij = y ( d j9u + diQji ~ d iOij) 

gjl 

y \pi9ij 4" digji diQij ) 

= 

= Y g 99 3l digji 
= —dig 

2 g 


(Eq. 105 with k = j) 

(relabeling dummy j, l in l si term & g = g /j ) 


(116) 


(derivative of determinant) 


In orthogonal coordinate systems, the Christoffel symbols of the hrst kind are given by: 


= DM] = o d j9v. 


no sum on i) 


[nJ] = ~2 d o9v 


(i 7^ j, no sum on i) (117) 

[ij, k] = 0 (i^jjLk) 

The hrst relation is a special case of Eq. 112 with l = i taking into account that the 
Christoffel symbols are symmetric in their paired indices; moreover, the relation includes 
the case of i — j, i.e. when all the three indices are identical. 

• In orthogonal coordinate systems, the Christoffel symbols of the second kind are given 
by: 


r* = 

1 jk 


DM] 


= g n DM] 


(no sum on i) 


(118) 
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and hence from the results of the previous point we have: 



(no sum on i) 



ii 



(no sum on i or j, and i 7^ j ) (H9) 



(»/j^ k) 


As for the first kind in the last point, the first relation includes the case of i — j , i.e. when 
all the three indices are identical. 

• In general orthogonal coordinate systems in a 3D space, the Christoffel symbols of the 
first kind vanish when the indices are all different, as shown earlier; moreover, the non- 
identically vanishing symbols of the first kind are given by: 


where (1,2,3) stand for (n 1 ,^ 2 ,^ 3 ) respectively, hi,h 2 ,h 3 are the scale factors as defined 
previously, and the comma indicates, as always, partial derivative; for example in cylin- 
drical coordinates given by (p,4>,z), h 2y 1 means the partial derivative of h 2 with respect 
to the first coordinate and hence h 2 ,\ = d p p = 1 since h 2 = p and the first coordinate is 
p (refer to Table 1). Because the Christoffel symbols of the first kind are symmetric in 
their first two indices, the [21, 1] symbol for instance can be obtained from the value of 
the [12, 1] symbol. 

• In general orthogonal coordinate systems in a 3D space, the Christoffel symbols of the 


[11.1] — -f/ii/ipi [11)2] — —hihip [11)3] — —h\h\£ 

[12. 1] = -\-h\h\p [12, 2] = T h 2 h 2 \ [13, 1] = - \-h\h\ ; 3 


[13, 3] = +h 3 h 3>1 [22,1 ] = -h 2 h 2>1 [22, 2] = +h 2 h 2)2 (120) 


[22,3] — —h 2 h 2 $ [23,2] — -\-h 2 h 2i 3 [23,3] — +h^h;^ 2 

[33, 1] = -h 3 h 3 , 1 [33, 2] = -h 3 h 3t 2 [33, 3] = +h 3 h 3 , 3 
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second kind vanish when the indices are all different, as shown earlier; moreover, the 
non-identically vanishing symbols of the second kind are given by: 


r 1 ■ 

h] A 

F 2 ■ 

/ll/ll, 2 

F 3 ■ 



h 

1 n - 

( h 2 ) 2 

1 li — 

(h 3 f 

r 1 — 

h 1,2 

h 

F 2 - 
1 12 - 

^2,1 
" " l_ h 2 

F 1 — 

1 13 — 

1-1 CO 

r 3 - 

h 3 ,i 

F 1 

h 2 h 2 ,i 

F 2 ■ 

^2,2 


+ h 

1 22 “ 

(hi ) 2 


h 2 

r 3 — 

^2^2,3 

0 hs ) 2 

F 2 - 

1 h 2,3 

' " l " h 2 

F 3 — 

h- 3,2 

h 3 

r 1 

h-ih-iA 

F 2 - 

/l 3 ^ 3,2 

F 3 ■ 

h 3j3 

1 33 " 

{hi? 

1 33 “ 

(h 2 f 

1 33 — 

+ h. 


( 121 ) 


where (1, 2, 3) stand for (m 1 ,^ 2 ,^ 3 ) respectively. Again, since the Christoffel symbols of the 
second kind are symmetric in their lower indices, the missing non-vanishing entries can be 
obtained from the given entries by permuting the lower indices. 

• In Cartesian coordinate systems ( x,y,z ), all the Christoffel symbols of the first and 
second kind are identically zero. 

• In cylindrical coordinate systems ( p,<t>,z ), the non- zero Christoffel symbols of the first 
kind are: 


[22.1] = -p (122) 

[ 12 . 2 ] = [ 21 , 2 ] = p 


where (1,2,3) stand for ( p,(j>,z ) respectively. 

• In cylindrical coordinate systems (p, (j), z ), the non-zero Christoffel symbols of the second 
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kind are: 


r 


1 

22 


-p 


r 


2 

12 



(123) 


where (1,2,3) stand for ( p,4>,z ) respectively. 

• In spherical coordinate systems (r, 9,(f>), the non- zero Christoffel symbols of the first 
kind are: 


[22,1] = 

— r 

[33,1] = 

— r sin 2 6 

[12,2] = 

[21,2]= r 

[33,2] = 

— r 2 sin 6 cos 6 

[13,3] = 

[31,3]= r sin 2 9 

[23,3] = 

[32, 3] = r 2 sin 9 cos 9 


(124) 


where (1,2,3) stand for ( r,9,(f > ) respectively. 

• In spherical coordinate systems (r,9,(f>), the non-zero Christoffel symbols of the second 
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kind are: 


r 

r 

r 

r 

r 

r 


1 

22 

— r 


1 

33 ~ 

—r sin 2 

9 

2 

12 ~ 

T 2 - 

1 

r 

CO to 
CO 

— sin 9 

cos 9 

3 

13 

r 3 — 

1 

r 

:i 

23 ~ 

r 3 — 

1 32 — 

cot i 


(125) 


where (1,2,3) stand for (r,9,(p) respectively. 

• Because there is an element of arbitrariness in the choice of the coordinates order and 
hence their indices, the Christoffel symbols may be given in terms of coordinate symbols 
rather than their indices to be more explicit and avoid ambiguity; for instance in the 
above examples of cylindrical and spherical coordinate systems we have: [22, 1] = [ff, p\, 
T 2 2 = rC for cylindrical, and [22, 1] = [99, r], T 3 3 = T^ for spherical. 

• The Christoffel symbols may be subscripted by the symbol of the metric tensor for the 
given space to reveal the metric which they are based upon. 

• In any coordinate system, all the Christoffel symbols of the first and second kind vanish 
identically iff all the components of the metric tensor in the given coordinate system are 
constants. 

• In affine coordinates, all the components of the metric tensor are constants and hence 
all the Christoffel symbols of both kinds vanish identically. 

• The number of independent Christoffel symbols of each kind (first and second) in general 
curvilinear coordinates is given by: 


N C i = 


n 2 (n + 1) 


2 


(126) 
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where n is the space dimension. The reason is that, due to the symmetry of the metric 
tensor there are rin+A independent metric components, and for each independent 
component there are n distinct Christoffel symbols. 

• The following relations are useful in the manipulation of tensor expressions involving 
Christoffel symbols: 16 


1 

77 


9i {V99 ij ) + g kl r« = o 


(127) 


dj [ik, l } = g la djV a ik + Y a ik [lj, a } + T“ fc [aj, l ] (128) 

3.2 Covariant Derivative 

• The basis vectors in general curvilinear coordinate systems undergo changes in magnitude 
and direction as they move around in their own space, and hence they are functions of 
position. These changes should be accounted for when calculating the derivatives of tensors 
in such general systems. Therefore, terms based on using Christoffel symbols are added 
to the ordinary derivative terms to correct for these changes and this more comprehensive 
form of derivative is called the covariant derivative. 

• Since in rectilinear coordinate systems the basis vectors are constants, the Christoffel 
symbol terms vanish identically and hence the covariant derivative reduces to the ordinary 
derivative, but in the other coordinate systems these terms are present in general. 

• As a consequence of the last point, the ordinary derivative of a non-scalar tensor is a 
tensor iff the coordinate transformations are linear. 

• It has been stated that the “covariant” label is an indication that the differentiation 
operator, V ; j, is in the covariant position. However, it may also be true that “covariant” 
means “invariant” as pointed out earlier in the previous set of notes. 

16 The first relation is a special case of the relation: A 1 ^ = ^jdj (yTh 4 lJ ) + A kl T kl noting that the 
covariant derivative of the metric tensor is identically zero according to the Ricci Theorem. 
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• Contravariant differentiation (V ;J ) can also be defined for covariant and contravariant 
tensors by raising the differentiation index using the index raising operator, e.g. 

A^=g^ k A i]k & A i »=g*A i . k (129) 

However, practically such operations are rarely used. 1 ' 

• As an example of how to obtain the covariant derivative of a tensor, let have a vector 
represented by contravariant components: A = A*E/ in general coordinates. We differ- 
entiate this vector following the normal rules of differentiation and taking account of the 
fact that the basis vectors in general coordinates are differentiable functions of position 
and hence they, unlike their Cartesian counterparts, are subject to differentiation using 
the product rule, that is: 

djA = E idjA 1 + A l djE, (product rule) 

= E .djA* + AT fc ..E, (Eq. 107) 

= E idjA 1 + A k Y l kj E, (relabeling dummy indices i & k) (130) 

= (djA* + A k r kj ) E* 

= 

where A* •, which is a rank-2 mixed tensor, is labeled the “covariant derivative” of A 1 . 
Similarly, for a vector represented by covariant components: A = A/E' in general curvi- 
linear coordinates we have: 

djA = Ai-jE* (131) 

• Following the method and techniques outlined in the previous point, to obtain the 
covariant derivative of a tensor in general, we start with an ordinary partial derivative 

1 ' An example of contravariant differentiation is in the definition of the Laplacian in general curvilinear 
coordinates (refer to § 4.1.4). 
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term of the given tensor. Then for each tensor index an extra Christoffel symbol term 
is added, positive for contravariant indices and negative for covariant indices, where the 
differentiation index is one of the lower indices in the Christoffel symbol. Hence, for a 
general differentiable rank-n tensor A the covariant derivative is given by: 


Aij...k 
^ lm...p\q 

f) A V- k 1 P* A a j-- k , -pi A ia - k 1 . , 
u <}- rs -lm...p ' 1 aq- rl lm...p ' 1 aq- tx lm...p ' 

<3 

£ — i 
+ 

(132) 


pa /pi...fc pa 

1 lq-^am...p 1 mq^la...p 

pa /pi--.fc 

L pq- n -lm...a 



• Practically, there is only one possibility for the arrangement of the indices in the Christof- 
fel symbol terms if the following rules are observed: 

(A) the second subscript index of the Christoffel symbol is the differentiation index, 

(B) the concerned tensor index in the Christoffel symbol term is contracted with one of 
the indices of the Christoffel symbol and hence they are opposite in their lower/upper 
position, 

(C) the contracted index is transferred from the tensor to the Christoffel symbol keeping 
its lower/upper position, and 

(D) all the other indices of the tensor keep their names and position. 

• The ordinary partial derivative term in the above covariant derivative expression (Eq. 
132) represents the rate of change of the tensor components with change of position as a 
result of moving along the coordinate curve of the differentiated index, while the Christoffel 
symbol terms represent the change experienced by the local basis vectors as a result of the 
same movement. This can be seen from the development of Eq. 130. 

• From the above discussion it is obvious that to obtain the covariant derivative, the 
Christoffel symbols are required and these symbols are dependent on the metric tensor; 
hence the covariant derivative is dependent on having the space metric. 

• In all coordinate systems, the covariant derivative of a differentiable scalar function of 
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position, /, is the same as the ordinary partial derivative, that is: 



( 133 ) 


This is justified by the fact that the covariant derivative is different from the ordinary 
partial derivative because the basis vectors in general coordinate systems are dependent 
on their spatial position, and since a scalar is independent of the basis vectors the covari- 
ant derivative and partial derivative are identical. This can also be concluded from the 
covariant derivative rule as stated in the previous points and formulated in Eq. 132. 

• Several rules of normal differentiation are naturally extended to covariant differentiation. 
For example, covariant differentiation is a linear operation with respect to algebraic sums 
of tensor terms and hence the covariant derivative of a sum is the sum of the covariant 
derivatives of the terms: 


where a and b are scalar constants and A and B are differentiable tensors. The product 
rule of ordinary differentiation also applies to covariant differentiation of inner and outer 
products of tensors: 


where the symbol o denotes an inner or outer product operator. 

• According to the “Ricci Theorem”, the covariant derivative of the covariant and con- 
travariant metric tensor is zero. This has nothing to do with the metric tensor being a 
constant function of coordinates, which is true only for the rectilinear systems, but this 
arises from the fact that the covariant derivative quantifies the change with position of the 
basis vectors in magnitude and direction as well as the change in components, and these 


(aA±f>B). j = a(A). j ±&(B) 


(134) 


(A o B) ;j = (A)., o B + A o (B) ;j 


(135) 
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contributions in the case of the metric tensor cancel each other resulting in a total null 
effect. As a result, the metric tensor behaves as a constant with respect to the covariant 
derivative operation: 

9ij-k = 0 & g.j. = 0 (136) 

for all values of the indices, and hence the covariant derivative operator bypasses the metric 
tensor: 

(g° A );fc = g° ( A );fc ( 13T ) 

where A is a general tensor, g is the metric tensor in its covariant or contravariant form 
and o denotes an inner or outer tensor product. 18 

• As a result of the Ricci Theorem, the covariant derivative operator and the index shifting 
operator are commutative, e.g. 


( 9 ikA k ).. = g lk A k - = Ai.j = (Ai).j = (< g ik A k ) ;j (138) 

( 9mi9njA mn ). k = gmi9njA = Aij- k = ( A ij). fe = ( g m ig n jA rnn ). k (139) 

• Like the metric tensor, the Kronecker delta is constant with regard to the covariant 
differentiation and hence the covariant derivative of the Kronecker delta is identically 
zero: 

S). k = d t s‘ + yry - s‘j% = o + ry _ t% = o (mo) 

The rule of the Kronecker delta may be regarded as an instance of the rule of the metric 
tensor, as stated by the Ricci Theorem, since the Kronecker delta is a metric tensor. 
Likewise, the covariant differentiation operator bypasses the Kronecker delta which is 

18 Although the metric tensor is normally used in inner product operations for raising and lowering of 
indices, the possibility of its involvement in outer product operations should not be ruled out. 
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involved in inner and outer tensor products: 19 


(<5oA) ;fc = <5o(A) ;fc (141) 

• Like the ordinary Kronecker delta, the covariant derivative of the generalized Kronecker 
delta is identically zero. 

• For a differentiable function f(x, y ) of class C 2 (i.e. all the second order partial derivatives 
of the function do exist and are continuous), the mixed partial derivatives are equal, that 
is: 

dxdyf = d y d x f (142) 

However, even if the components of a tensor satisfy this condition (i.e. being of class C 2 ), 
this is not sufficient for the equality of the mixed covariant derivatives. What is required 
for the mixed covariant derivatives to be equal is the vanishing of the Riemann Tensor 
(see § 5.1). 

• Higher order covariant derivatives are similarly defined as derivatives of derivatives 
by successive repetition of the process of covariant differentiation; however the order of 
differentiation should be respected as stated in the previous point. For example, the second 
order mixed jk covariant derivative of a contravariant vector A is given by: 


A'vk = dkdjA' 


+ n a d 3 A° 


jk 


d a A l + T l ja d k A c 


+ A a (d 


r* 

J L ka 


r b r* -i- r 

1 jk L k + i 


jXL) (143) 


while the second order mixed jk covariant derivative of a covariant vector A is given by: 


Ai-jk dkdjAi 


r ^d k A a 




r - A (Ar“ 


r a r 


ib L jk 


-pb -pa \ 

ik^~ bj ) 


(144) 


19 Like the metric tensor, the Kronecker delta is normally used in inner product operations for replace- 
ment of indices; however the possibility of its involvement in outer product operations should not be ruled 
out. 
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• The second order mixed kj covariant derivative of contravariant and covariant vectors 
can be obtained from the equations in the last point by interchanging the j and k indices 
and hence the inequality of the jk and kj mixed derivatives in general can be verified 
(refer to § 5.1). 

• The covariant derivative of a tensor is a tensor whose covariant rank is higher than the 
covariant rank of the original tensor by one. Hence, the covariant derivative of a rank-n 
tensor of type (r, s ) is a rank-(n + 1) tensor of type (r, s + 1). 

• Covariant differentiation and contraction of index operations commute with each other, 
e.g. 

(4VU(4U ; , = (4*) ; , ( 145 > 

• Since the Christoffel symbols vanish when the components of the metric tensor in a 
given coordinate system are constants, the covariant derivative is reduced to the ordinary 
derivative in such systems. This is particularly true in Euclidean spaces coordinated by 
rectilinear systems. 

• The covariant derivatives of relative tensors, which are also relative tensors of the same 
weight as the original tensors, are obtained by adding a weight term to the normal formulae 
of covariant derivative. Hence, the covariant derivative of a relative scalar with weight w 
is given by: 

./': > (146) 

while the covariant derivative of relative tensors of higher ranks with weight w is obtained 
by adding the following term to the right hand side of Eq. 132: 

-«4kb4 < 147 > 

• Unlike ordinary differentiation, the covariant derivative of a non-scalar tensor with con- 
stant components is not zero in general due to the presence of the Christoffel symbols in 
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the definition of the covariant derivative, as given by Eq. 132. 

• In rectilinear coordinates, the Christoffel symbols are identically zero because the basis 
vectors are constants, and hence the covariant derivative is the same as the normal partial 
derivative for all tensor ranks. As a result, when the components of the metric tensor, , 
are constants as in the case of rectangular coordinate systems, the covariant derivative 
becomes the ordinary partial derivative. 

• For a differentiable covariant vector A which is a gradient of a scalar we have: 


• The covariant derivative of the basis vectors of the covariant and contravariant types is 
identically zero: 


3.3 Absolute Derivative 

• The absolute derivative of a tensor along a t-parameterized curve in an nD space with 
respect to the parameter t is the inner product of the covariant derivative of the tensor and 
the tangent vector to the curve. In brief, the absolute derivative is a covariant derivative 
of a tensor along a curve. 

• For a tensor A\ the inner product of A 1 -, which is a tensor, with another tensor is a 
tensor. Now, if the other tensor is which is the tangent vector to a f-parameterized 
curve L given by the equations u l = u l (t ), then the inner product: 


A . . — A . . 
3 A?;* 


(148) 


F,, = djEj - I 'f F- = Tj-Efc - r'; ; F/, = 0 
Eh = <),e : + r kj E k = -r fcj E fe + r;, e ,v = o 


(149) 



(150) 


is a tensor of the same rank and type as the tensor Ah The tensor given by the expression 
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150 is called the “absolute” or “intrinsic” or “absolute covariant” derivative of the tensor 
A 1 along the curve L and is symbolized by: 


8 A 

IT 


(151) 


• For a differentiable scalar /, the absolute derivative, like the covariant derivative, is the 
same as the ordinary derivative, that is: 


b_i = df_ 

5t dt 


(152) 


• The absolute derivative of a differentiable contravariant vector A k with respect to the 
parameter t is given by: 


5A^ = dA k fc jdvP 

5t dt y dt 


(153) 


Similarly for a differentiable covariant vector Ak we have: 


5A k _ dA k ■ du 3 

5t ~ dt kj 1 dt 


(154) 


• Absolute differentiation can be easily extended to higher rank (> 1) differentiable tensors 
of type (m, n) along parameterized curves. For instance, the absolute derivative of a mixed 
tensor of type (1,2) A- k along a /-parameterized curve L is given by: 


St 



du b 

dt 


dA 


jk 


dt 


+ r 


ab 


A a ^ 

jk dt 


F“ A — 
jb ak dt 


F« 4* d A 

hk ja dt 


(155) 


• As the absolute derivative is given generically by: 


= (A) ;fc 


du k 

dt 


6 A 

~6t 


(156) 
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it can be seen as an instance of the chain rule of differentiation where the two contracted 
indices represent the in-between coordinate differential. 

• Because the absolute derivative along a curve is just an inner product of the covariant 
derivative with the vector tangent to the curve, the well known rules of ordinary dif- 
ferentiation of sums and products also apply to absolute differentiation, as for covariant 
differentiation, that is: 

i( tt A + ®) = «f (157) 


!< A ° B )=(f ° B M A °f) (158) 

where a and b are constant scalars, A and B are differentiable tensors and the symbol o 
denotes an inner or outer product of tensors. 

• The covariant and contravariant metric tensors are in lieu of constants with respect to 
absolute differentiation, that is: 


SQij 

St 


0 


& 



(159) 


and hence they pass through the absolute derivative operator: 


S (gijA 3 ) _ 5AJ 

— gij- 


st 


st 


& 


Hg ij Aj 

St 


= 9 if 


(160) 


• For coordinate systems in which all the components of the metric tensor are constants, 
the absolute derivative is the same as the ordinary derivative. This is the case in the 
rectilinear coordinate systems. 

• The absolute derivative of a tensor along a given curve is unique, and hence the ordinary 
derivative of the tensor along that curve in a rectangular coordinate system is the same 
as the absolute derivative of the tensor along that curve in any other system. 
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4 Differential Operations 

• In this section we generalize and expand what have been given in the previous notes 
[11] about the main differential operations which are based on the nabla operator V. The 
section will investigate these operations in general curvilinear coordinate systems and in 
general orthogonal coordinate systems which are a special case of the general curvilinear 
systems. We also investigate the two most important and widely used non-Cartesian 
orthogonal coordinate systems, namely the cylindrical and spherical systems, due to their 
particular importance. 

4.1 General Curvilinear Coordinate System 

• Here, we investigate the differential operations and operators in general curvilinear 
coordinate systems, whether orthogonal or not. 

• The previous definitions of the differential operations, as given in the first set of notes, 
are essentially valid in general non-Cartesian coordinate systems if the operations are 
extended to include the basis vectors as well as the components. 

• The analytical expressions of the differential operations can be obtained directly if the 
expression for the nabla operator V and the spatial derivatives of the basis vectors in the 
general curvilinear coordinate system are known. 

4.1.1 Gradient 

• The nabla operator V in general curvilinear coordinate systems is defined as follow: 


V = E'di 


(161) 


4.1.1 Gradient 


60 


Hence, the gradient of a differentiable scalar function of position, /, is given by: 


V/ = E 'dif = E'/i 


(162) 


The components of this expression represent the covariant form of a rank-1 tensor, i.e. 
[V/] . = j \ t , as it should be since the gradient operation increases the covariant rank 
of a tensor by one. Since this expression consists of a contravariant basis vector and 
a covariant component, the gradient in general curvilinear systems is invariant under 
admissible transformations of coordinates. 

• The contravariant form of the gradient of a scalar / can be obtained by using the index 
raising operator, that is: 


• The gradient of a differentiable covariant vector A can similarly be defined as follow: 


[ v /] 1 = <97 = gVdjf = g'-'.f.i = f 


(163) 


VA = E*<9j (AjE j ) 


= E ‘E’diAj + E i A j d i E j 


(product rule) 


= E E (9 ; A j + E i Aj (-T j ki E k ) 


(Eq. 108) 


(164) 


= E'E J di A j - E'E'Y^Ak 
= E Iv (diAj - Y%A k ) 

= E E Aj.j 


(relabeling dummy indices j & k) 


(taking common factor E'E ; ) 


(definition of covariant derivative) 
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Similarly, for a differentiable contravariant vector A the gradient is given by: 


VA = E % {A j Ej) 

= E'E jdiA j + E'.V(7;E: 

= E'E ,diA j + EM' ( I 'j; E/, ) 
= E'E ,diA j + E%Y j ki A k 
= E-E, (a.v + 4^) 

= E'E,/y 


(product rule) 

(Eq. 107) 

(165) 

(relabeling dummy indices j & k ) 

(taking common factor E'E;) 

(definition of covariant derivative) 


The components of the gradients of covariant and contravariant vectors represent, respec- 
tively, the covariant and mixed forms of a rank-2 tensor, as they should be since the 
gradient operation increases the covariant rank of a tensor by one. 

• The gradient of higher rank tensors is similarly defined. For example, the gradient of a 
rank-2 tensor is given by: 


VA 

= E'E'E" (diA jk - Y l 3i A lk - T l ki Aji) = E'E'E k A jk , i 

(covariant) 

(166) 

VA 

= E , E J E k ( dtA jk + Y j u A lk + Y k u A jl ) = E%E k A jk 

(contravariant) (167) 

VA 

= E'E'E, (diA k - T^Af + t£A/) = E'E'E.A^ 

(mixed) 

(168) 

VA 

= E*EjE fc (d t A\ + Y] t A l k - Y l kl A \ ) = E i EjE k A\. i 

(mixed) 

(169) 


4.1.2 Divergence 

• Generically, the divergence of a differentiable contravariant vector A is defined as follow: 


V ■ A = E l di ■ (A j Ej) = E' ■ di (A%) 


e‘ • (y %) = (E‘ ■ e 3 ) a% = ajy, 


A% (170) 


4.1.2 Divergence 


62 


In more details, the divergence of a differentiable contravariant vector A 1 is a scalar ob- 
tained by contracting the covariant derivative index with the contravariant index of the 
vector, and hence: 


V • A = A\; 


= diA 1 + t;,tb 

= d i A i + Ai^-d j { y /g) 
V9 

= d . A i + A iJ_ d .^) 

Vg 


(definition of covariant derivative) 
(Eq. 115) 

(renaming dummy index j) 


(171) 


(product rule) 


where g is the determinant of the covariant metric tensor g %3 . The last equality may be 
called the Voss-Weyl formula. 

• The divergence can also be obtained by raising the first index of the covariant derivative 
of a covariant vector using a contracting contravariant metric tensor: 


g^Aj;, = (g s 'A,). f 

(Eq. 137) 

= (A% 

(raising index) 

= A‘. i 


= V-A 

(Eq. 171) 


(172) 


as before. 

• Based on the previous point, the divergence of a covariant vector Aj is obtained by using 
the raising operator, that is 

A], = g^Aj-i 


(173) 
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• For a rank- 2 contravariant tensor A, the divergence is generically defined by: 


V • A = E% • (A j %E k ) = (E • E y ) E k A% = 8 i j E k A j ^ i = E k A% (174) 

The components of this expression represent a contravariant vector, as it should be since 
the divergence operation reduces the contravariant rank of a tensor by one. 

• More generally, considering the tensor components, the divergence of a differentiable 
rank- 2 contravariant tensor A l i is a contravariant vector obtained by contracting the co- 
variant derivative index with one of the contravariant indices, e.g. 

[V • A] j = A ij - (175) 

And for a rank- 2 mixed tensor A 1 - we have: 

[V . A], = 4, (176) 

• Similarly, for a general tensor of type (m,n): A = the divergence with 

respect to its k th contravariant index is defined by: 

[v-a]‘s:::4 = (4|v:::4) ; , ( 177) 

with the absence of the contracted contravariant index i k on the left hand side. 
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4.1.3 Curl 

• The curl of a differentiable vector is the cross product of the nabla operator V with the 
vector, e.g. the curl of a vector A represented by covariant components is given by: 


curl A = V x A 


= E % x AjEP 
= E* x di (Aj E j ) 
= E* x (Aj.0) 
= A j;i (E* x E- 7 ) 
= Aj-i e ijk E k 


(Eq. 131) 


(Eq. 78) 


(178) 


= e ijk A j;i E h 

ijk 

= (Mj - rV,) e„ 


(Eq. 49 & definition of covariant derivative) 


and hence the k th contravariant component of curl A is: 


. Jjk 

[V x A]* = — (dtAf - r CA) (179) 

On expanding the last equation for the three components of a 3D space, considering that 
the terms of the Christoffel symbols cancel out due to their symmetry in the two lower 
indices, 20 we obtain: 


[V x A] 1 
[V x A] 2 
[V x A] 3 


((? 2 x 43 — d^A^) 

V9 

—p (d^Ai — d\A$) 

V9 

—— {d\A<2 — dzAi) 

V9 


20 


That is: — Aj — djA^ — A^T^ — d{Aj + A^T^ — djA 1 — A^T^j — diAj + A^T^a 


A ■ _ A ■ ■ 


(180) 

(181) 

( 182 ) 
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Hence, Eq. 179 will reduce to: 


[V x Af 


ijk 

-p diAj 

V9 


(183) 


4.1.4 Laplacian 

• Generically, the Laplacian of a differentiable scalar function of position, /, is defined as 
follow: 

V 2 / = div (grad /) = V • ( V/) (184) 


Hence the simplest approach for obtaining the Laplacian in general coordinates is to insert 
the expression for the gradient, V/, into the expression for the divergence. However, 
because in general curvilinear coordinates the divergence is defined only for contravariant 
tensors whereas the gradient of a scalar is a covariant tensor, the index of the gradient 
should be raised first before applying the divergence operation, that is: 


[v/]‘ = 57 = 9 «a ,7 


(185) 


Now, according to Eq. 171 the divergence is given by: 


V ■ A = A i . i 


1 

V5 


Si (s/gA‘) 


(186) 


On defining A = V/ 


E id 1 f and replacing A 1 in the last equation with d l f we obtain: 


V7 = V • (E,57) = —pSi (VggSdjf) (187) 

V9 

which is the expression for the Laplacian of a scalar function / in general curvilinear 
coordinates. 

• Another approach for developing the Laplacian expression in general coordinates is to 
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apply the first principles by using the definitions and basic properties of the operations 
involved, that is: 

v 2 / = V-(V/) 

= E% ■ (Wdjf) 

= i : ■ d; (E'djf) 

= E i -d i (Wf J ) 

= E' • (E Pfjf) 

= (E ; • E J ) fa 
= iTfM 
= (9V,i), 

= GW), 


d< (a%f) + (a^f) it. 

(Eq. 132) 

di (s%f) + (9%/) rt 

(renaming dummy indices i & k) 

9 . (yy/) + -L (a,v9) 

(Eq. 115) 

[Vsa (a%f) + (g%f) d iVa] 

(taking factor) 

y/9 

(Va^Sif) 

(product rule) 


as before. 

• The Laplacian of a scalar / may also be shorthand notated with: 

V 2 / = (189) 


(Eq. 131) 


(Eq. 51) 
(Eq. 137) 


(188) 


• The Laplacian of non-scalar tensors can be similarly defined. For example, the Laplacian 
of a vector B is a vector A (i.e. A = V 2 B) which may be defined in general coordinates 
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as: 

A* = and Ai = g jk B i]jk (190) 

• The Laplacian of a tensor is a tensor of the same rank and variance type. 


4.2 General Orthogonal Coordinate System 


• In this section we state the main differential operations in general orthogonal coordinate 
systems. These operations are special cases of the operations in general curvilinear systems 
which were derived in the previous section. However, due to the wide spread use of 
orthogonal systems, it is worth to state the most important of these operations although 
they can be easily obtained from the formulae of general curvilinear systems. 

• General orthogonal coordinate systems are identified in the following notes by the co- 
ordinates (u 1 , . . . , u n ) with unit basis vectors (u 1; . . . , u„) and scale factors (hi, . . . , h n ) 
where: 21 

v - -' 1 Ox 3 \ — ' 0v7 . . . . 

Uj = > — — -e,- = > m— — re,- (no sum on i) (191) 

^ hidu 1 3 ^ Ox 3 3 v ’ v ’ 

3 3 


K = IE; I = 


Or 


Out 


f Ox 3 \ ' 

\ du l J 


L 3 


1/2 


E 


du l 

Ox 3 


- 1/2 


(192) 


In the last equations, x 3 and e 7 - are respectively the coordinates and unit basis vectors in 
the Cartesian rectangular system, and r is the position vector in that system. 


4.2.1 Gradient 

• The nabla operator in general orthogonal coordinates is given by: 


v = E 

i 


U; 0 

hi dv7 


(193) 


21 In orthogonal coordinates, the covariant and contravariant normalized basis vectors are identical, as 
established previously in § 1.7, and hence u* = u, ; and e 3 = ej. 
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Hence, the gradient of a differentiable scalar / in orthogonal coordinates is given by: 


v/ = £ 

i 


lb df 
hi du i 


(194) 


4.2.2 Divergence 


• The divergence in orthogonal coordinates can be obtained from Eq. 171. Since for 
orthogonal coordinate systems the metric tensor is diagonal with ^fg = h\h 2 h 3 in a 3D 
space and hiA 1 = A 1 (component-wise with no summation), the last line of Eq. 171 
becomes: 


V- A 


i_d_ 

Jgdu 


(v^) 


1 d / h\h 2 h 3 

hih 2 h 3 du i \ hi 


(195) 


where A is a contravariant differentiable vector and A2 represents the physical components 
(refer to Eq. 28). 22 This equation is the divergence of a vector in general orthogonal 
coordinates as defined in vector calculus. 


4.2.3 Curl 


• The curl of a differentiable vector A in orthogonal coordinate systems in 3D spaces is 
given by: 

hi Ui h 2 u 2 h 3 u 3 


V x A 


1 

hih 2 h 3 


d 
du 1 


a a 

du 2 du 3 


(196) 


h\Ai h 2 A 2 h 3 A 3 


where the hat indicates a physical component. The last equation may also be given in a 
more compact form as: 


[V7 ^ A1 _ \ " e ijkhj d(h k A k ) 
^ hih 2 h 3 dui 

k = 1 


(197) 


22 In orthogonal coordinate systems the physical components are the same for covariant and contravari- 
ant forms, as established before in § 1.7, and hence Aj = A*. 
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4.2.4 Laplacian 


• For general orthogonal coordinate systems in 3D spaces we have: 


y/g = hih 2 h 3 & g lJ = — — (no sum on i or j) 

I Li I In 


(198) 


and hence Eq. 187 becomes: 



(199) 


which is the Laplacian of a scalar function of position, /, in orthogonal coordinates as 


defined in vector calculus. 

4.3 Cylindrical Coordinate System 

• For cylindrical systems identified by the coordinates (p,(j),z), the orthonormal basis 
vectors are e p , and e,. 23 We use for brevity e pr p as a shorthand notation for the unit 
dyad and similar notations for the other dyads. 

4.3.1 Gradient 

• The gradient of a differentiable scalar / is: 


23 Hence the given components (i.e. and A z ) are physical (see § 1.7). Despite that, we do not 

use hats since the components are suffixed with coordinate symbols (refer to § 1.7). 



( 200 ) 
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The gradient of a differentiable vector A is: 


VA GppApp T g ptyA^p T GpzA z ^p T 

*4 >p 

^zpAp^z T T G zzA Z: z 


. 1 . 7L\ (1 A A p \ 1 4 

e <t>o I -A p < f, + e^ —A^ fj, H + e^z—Az t/, + 

P P ) \P P ) P 


(201) 


4.3.2 Divergence 


The divergence of a differentiable vector A is: 


V • A — - [dp ( pA p ) + d^A^ + pd z A z ] 


(202) 


The divergence of a differentiable rank-2 tensor A is a vector given by: 


24 


V-A = 


e 0 ( A PPt p + App ^ Ah + l -A^ + A zp ,z ) + 


I A-p<t>,p T 


2A 


p<I> 


H — + Azfc- + 


P P 


■At j>p A P0 

P 


e z ( A P z tP H — + -A^ z ,(j> + A zz , z 


+ 


(203) 


4.3.3 Curl 


The curl of a differentiable vector A is: 


V x A = - 
P 


® p P e (f) 


dp d# d z 


Ap pA^ A z 


(204) 


24 It should be understood that V ■ A is lower than the original tensor A by just one contravariant 
index and hence, unlike the common use of this notation, it is not scalar in general. 
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4.3.4 Laplacian 

• The Laplacian of a differentiable scalar / is: 


V 2 / = d pp f + -d p f + + d zz f 

p /r 


(205) 


The Laplacian of a differentiable vector A is: 


V 2 A = e £ 


dp (^—®p + -^d^Ap + d zz A p — —d^A^ 

dp (^—^p + -^d^Afj, + d^A^ + —d^Ap 

—dp ( pd p A z ) + —dtj^Az + d zz A z 


+ 

+ 


(206) 


4.4 Spherical Coordinate System 

• For spherical coordinate systems identified by the coordinates (r,9,<p), the orthonormal 
basis vectors are e r , e g and e^. 25 We use for brevity e r9 as a shorthand notation for the 
dyad e r eo and similar notations for the other unit dyads. 


4.4.1 Gradient 

• The gradient of a differentiable scalar / is: 

V/ = e r d r f + e g -d g f + <9 0 / 

r r sm t) 


(207) 


25 Again, the components are physical and we do not use hats. 
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The gradient of a differentiable vector A is: 


VA ^^yy y A G r eAo t r A ^r(f>A(p^ r A 


( 208 ) 


Gg r 


-•(jyr 


A, 


r,6 


r 

A, 


r,<f> 


r sin 9 


A e 

— ) + e e9 
r 

A <), 

) + e d 

r 


A 9 a A r 


H I + ee<p- 


A 


(f>,8 


A 


A e ^ A^ cot 9 


r sin 9 


+ e,^ 


A^ ^ A r ^ Aq cot 6 


r sin 9 r 


4.4.2 Divergence 


The divergence of a differentiable vector A is: 


V- A = 


1 


r 2 sin 9 


d ( r 2 A r ) d (sin 9A e ) d A, 

sin 9 „ ’ + r — - — — — — + r 


dr 


89 


d<p 


(209) 


The divergence of a differentiable rank-2 tensor A is a vector given by: 


V-A = 


d r (r 2 A rr ) d 9 (Ae r sin 9) d^A^ A ee + A H 
e, I o 1 ^ 1 t 77 + 


r sin 9 


r sin 9 


( 210 ) 


' d r (r 3 A rd ) d e (A ee sin 6) d^A^e A 0r - A r6 - A H cot 9 \ 
e e I ; 1 — 1 — H I + 


r sin 9 


r sin 9 


8 r (r A r( p ) 8q (Aecp sind) d^A^^ A<p r A r( p A A^q cot 9 

e rh I 5 1 : — 77 1 : — 77 H 


r sin 9 


r sin 9 


4.4.3 Curl 


The curl of a differentiable vector A is: 


V x A = 


r 2 sin 9 


e r re g r sin de^, 

d r do d 0 

A r rAe r sin 9A 4 


( 21 !) 
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5 Tensors in Applications 

• In this section we conduct a preliminary investigation of some commonly-used tensors in 
physical and mathematical applications of tensor calculus in anticipation of the forthcom- 
ing set of notes. Most of these tensors come from differential geometry, fluid, continuum 
and relativistic mechanics, since these disciplines are intimately linked to tensor calculus as 
large parts of the subject were developed, and are still developing, within these disciplines. 
These tensors also form the building blocks of several physical and mathematical theories. 
We would like to insist that these are just a few partially representative examples for the 
use of tensors in scientific and mathematical applications to have more familiarity with 
tensor language and techniques and hence they are not meant to provide a comprehensive 
view in any way. We should also indicate that some tensors are defined differently in 
different disciplines and hence the given definitions and properties may not be thorough 
or general. 

5.1 Riemann Tensor 

• This rank-4 tensor, which is also called Riemann curvature tensor and Riemann-Christoffel 
tensor, is a property of the space. It characterizes important properties of spaces and sur- 
faces and hence it plays an important role in geometry in general and in non-Euclidean 
geometries in particular. 

• The covariant differentiation operators in mixed derivatives are not commutative and 
hence for a covariant vector A we have: 

Aj-m Oj-ij, Rj k iAi 


(214) 
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where Kj kl is the Riemann tensor of the second kind which is given by: 



o k n 


jl 


r)V l I V r V 1 V r V 1 
jk jl*- rk jk^ rl 


(215) 


• The last equation can be put into the following mnemonic determinantal form: 



d k 

di 

+ 

f-i 

pr 

1 jk 

-pi 

1 jk 

L ji 


r* 

1 rl 

r* 

1 rk 


(216) 


• The Ricmann-Christoffel tensor of the second kind is also called the mixed Riemann- 
Christoffel tensor. 

• From Eq. 214, it is obvious that the mixed second order covariant derivatives are equal 
iff the Riemann tensor of the second kind vanishes identically. 

• On lowering the contravariant index of the Riemann tensor of the second kind, the 
Riemann tensor of the first kind is obtained: 


Rijkl — QiaR jkl 


(217) 


• Alternatively, the Riemann tensor of the first kind can be established independently as: 


Rijki = d k [jl, i] - di [jk, i] + [il, r] T r jk - [ ik , r] TJ, 

= \ ( djd k gn + didig jk - djdig ik - did k gji ) + [il, r] T r jk - [ik, r] ( 218 ) 

= ^ ( djdkgn + didig jk - djdig ik - d^gji) + g rs ([il, r } [jk, s } - [ik, r] [jl, s]) 


• The first line of the last equation can be cast in the following mnemonic determinantal 
form: 


R 


ijkl 


d k 

di 

+ 

pr pr 

1 jk 1 jl 

b'M 

WA 


[■ ik,r\ [il,r\ 


(219) 
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• Similarly, the Riemann-Christoffel tensor of the second kind can be obtained by raising 
the first covariant index of the Riemann-Christoffel tensor of the first kind: 

rtju = S ,a R*j a (220) 

• The Riemann-Christoffel tensor of the first kind is also called the covariant (or totally 
covariant) Riemann-Christoffel tensor. 

• The covariant differentiation operators become commutative when the metric makes the 
Riemann tensor of either kind vanish. 

• For the mixed second order covariant derivatives of a contravariant vector A we have: 

A\ki ~ A j - } ik — R?ukA l ( 221 ) 

which is similar to Eq. 214 for a covariant vector A. 

• The Riemann-Christoffel tensor vanishes identically iff the space is globally flat. Hence, 
the Riemann tensor is zero in Euclidean spaces, and consequently the mixed second order 
covariant derivatives, which become ordinary derivatives, are equal when the C 2 continuity 
condition is satisfied. 

• The Riemann curvature tensor depends only on the metric which, in general curvilinear 
coordinates, is a function of position and hence the Riemann tensor follows this dependency 
on position. Yes, for affine coordinates the metric tensor is constant and hence the Riemann 
tensor vanishes identically. 

• The totally covariant Riemann tensor satisfies the following symmetric and skew-symmetric 


5.1 Riemann Tensor 


77 


relations in its four indices: 


Rijki = Rklij (block symmetry) 

= —Rjiki (anti-symmetry in the first two indices) 

= — Rijik (anti-symmetry in the last two indices) 


( 222 ) 


• The skew-symmetric property of the covariant Riemann tensor with respect to the last 
two indices also applies to the mixed Riemann tensor: 


tdi m 

- a jkl — jlk 


(223) 


• As a consequence of the first and second anti-symmetric properties of the covariant 
Riemann tensor, the entries of the Riemann tensor with identical values of the first two 
indices or/and the last two indices are zero. 

• As a consequence of the skew-symmetric properties of the Riemann tensor, all entries of 
the tensor with identical values of more than two indices (e.g. Ruji) are zero. 

• In an nD space, the Riemann tensor has n 4 components. 

• As a consequence of the symmetric and anti-symmetric properties of the Riemann tensor, 
in an nD space there are three types of distinct non-vanishing entries: 

A. Entries with only two distinct indices (type Rijij) which count: 


Ah = 


n (n — 1) 


(224) 


B. Entries with only three distinct indices (type Rijik) which count: 


Ah = 


n (■ n — 1) (n — 2) 


(225) 
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C. Entries with four distinct indices (type R. l3 ki) which count 

Y n ( n ~ 1) ( n ~ 2) { n ~ 3) 

3 “ 12 

• By adding the numbers of the three types of non-zero distinct entries, as given in the 
last point, it can be shown that the Riemann tensor in an nD space has a total of 

rR (rR — R 

N m = N 1 + N 2 + N 3 = — ^ >- (227) 

independent components which do not vanish identically. For example, in a 2D Ric- 
mannian space the Riemann tensor has 2 4 = 16 components; however there is only one 
independent component (with the principal suffix 1212) which is not identically zero rep- 
resented by the following four dependent components: 

R\ 2 i 2 — R 2121 — ~R \221 — ~ -R 2112 (228) 

Similarly, in a 3D Riemannian space the Riemann tensor has 3 4 = 81 components but only 
six of these are distinct non-zero entries which are the ones with the following principal 
suffixes: 

1212, 1313, 1213, 2123, 3132, 2323 (229) 

where the permutations of indices in each of these suffixes are subject to the symmetric 
and anti-symmetric properties of the four indices of the Riemann tensor, as in the case of a 
2D space in the above example, and hence these permutations do not produce independent 
entries. 

• Following the pattern in the last point, in a 4D Riemannian space the Riemann tensor 
has 4 4 = 256 components but only 20 of these are independent non-zero entries, while in 
a 5D Riemannian space the Riemann tensor has 5 4 = 625 components but only 50 are 


(226) 
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independent non-zero entries. 

• The Riemann tensor satisfies the following identity: 

R-I]kl:s T Riljk-,s iksl;j T Rikjs\l (230) 

• A necessary and sufficient condition that a manifold for which there is a coordinate 
system with all the components of the metric tensor being constants 26 is that: 


Rij kl 0 


(231) 


• On contracting the first covariant index with the contravariant index of the Riemann 
tensor of the second kind we obtain: 


pi pi pi pr pi pr pi 

ikl VkL n C//1 ifc ~T L j^L 1 ifcl r i 

d pi pi I pr pi pi pr 

k*~ il ”1*- jjj il r /j. 1 n 

= d - d,r ik 


= dk [di (111 y/g)] - di [d k (In y/g)] 
= d k d i (In y/g) - d{d k (In y/g) 

= d k d i (In y/g) - d k di (In y/g) 


= 0 


(j = i in Eq. 215) 

(relabeling dummy i,r in last term) 

(Eq. 115) (232) 

(C 2 condition is assumed) 


That is: 

R'iki = o 


( 233 ) 


26 This may be called flat pseudo Riemannian manifold. 
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5.1.1 Bianchi identities 

• The Riemann tensor of the first and second kind satisfies a number of identities called 
the Bianchi identities. 

• The first Bianchi identity is: 


Rijki + Riijk + Rikij = 0 (first kind) 

R l jki + R\jk + R-'ki] — 0 (second kind) 


(234) 


These two forms of the first identity can be obtained from each other by the raising and 
lowering operators. 

• The above first Bianchi identity is an instance of the fact that by fixing the position of 
one of the four indices and permuting the other three indices cyclically, the algebraic sum 
of these three permuting forms is zero, that is: 


Rijki T Rujk T Rikij 0 ( i fixed) 

Rijki T Rljik + Rkjli 0 ( j fixed) 

Rijki ~fi Rlikj Rjlki d ( 1' fixed) 

Rijki T Rkijl + Rjkil 0 ( i fixed) 


• Another one of the Bianchi identities is: 

Rijki ;m + Rijhn\k T R'ijm,k:i 0 (first kind) 

R'dki-.m + R l jim,k + R’jmk/.i = 0 (second kind) 


(235) 


(236) 


Again, these two forms can be obtained from each other by the raising and lowering 
operators. 

• The Bianchi identities are valid regardless of the metric. 
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5.2 Ricci Tensor 

• The Ricci tensor of the first kind is obtained by contracting the contravariant index with 
the last covariant index of the Ricmann tensor of the second kind, that is: 

Rn = Rija = SjTi - d.r% + ny*. - rj„rt, (237) 


and hence it is a rank- 2 tensor. 

• The Ricci tensor, as given by the last equation, can be written in the following mnemonic 
determinantal form: 


% 

d a 

+ 

r a 

1 bj 

-na 

ba 

T“; 

T“ 


rk 

r b n 

*7 

ia 


0 

ia 


(238) 


• Because of Eq. 115 (i.e. = <9j (ln^)), the Ricci tensor can also be written in the 

following forms as well as several other forms: 


Rij = djdi (In y/g) - d a T “■ + T a bj T b a - T^db (In y/g) 

= dA (In y/g) + T a bj T b ia - d a T “ - T ^ (In y/g) (239) 

= djdi (In y/g) + T^T b a - d a (y/gT^) 

where g is the determinant of the covariant metric tensor. 

• The Ricci tensor of the first kind is symmetric, that is: 


Rij = Rji (240) 

• On raising the first index of the Ricci tensor of the first kind, the Ricci tensor of the 
second kind is obtained: 

R‘i = g ik Rkj 


(241) 
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• The Ricci scalar, which is also called the curvature scalar and the curvature invariant, 
is the result of contracting the indices of the Ricci tensor of the second kind, that is: 

R = R\ (242) 


• Since the Ricci scalar is obtained by raising a subscript index of the Ricci tensor of 
the first kind using the raising operator followed by contracting the two indices, it can be 
written as: 


R = g"R ; , = g lJ 


dA (In y/g) + T a bj T\ a - —d a (^ ) 


(243) 


where the expression in the square brackets is obtained from the last line of Eq. 239; 
similar expressions can be obtained from the other lines of that equation. 

• In an nD space, the Ricci tensor has n 2 entries. However, because of its symmetry it 
possesses a maximum of 


Ard — 


n (n + 1) 


(244) 


distinct entries. As an example, in the 4D manifold of general relativity n = 4, and hence 
the Ricci tensor has 4 2 = 16 components. However, due to the symmetry of the Ricci 
tensor there are only ten independent entries according to Eq. 244. The gravitational 
field equations in a free space are obtained by setting the Ricci tensor components equal 
to zero, and hence there are ten partial differential equations describing the gravitational 
field in this space according to the general relativistic mechanics. 
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5.3 Einstein Tensor 


• The Einstein tensor G is a rank-2 tensor defined in terms of the Ricci tensor R and the 


Ricci curvature scalar R as follow: 2 ' 


C 


c 



(covariant) 


(mixed) 


(contravariant) 


(245) 


• Since both the Ricci tensor and the metric tensor are symmetric, the Einstein tensor is 
symmetric as well. 

• The divergence of the Einstein tensor vanishes at all points of the space for any Rieman- 
nian metric. 

• On contracting the Bianchi identity twice with using the anti-symmetric properties of 
the Riemann tensor we obtain: 


which is inline with the above statement. The following form can also be derived based 
on the Bianchi identity: 


G m ™ = 0 


t mn 


(246) 


g jk G ki .j = 0 


(247) 


5.4 Infinitesimal Strain Tensor 


• This is a rank-2 tensor which describes the state of strain in a continuum medium and 


hence it is used in continuum and fluid mechanics. 


27 We notate this tensor with G rather than E, which is more natural, to avoid potential confusion with 
the first displacement gradient tensor (see § 5.6). 
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• The infinitesimal strain tensor 7 is defined by : 28 


7 


Vd + Vd T 

2 


(248) 


where d is the displacement vector and the superscript T represents matrix transposi- 
tion. The displacement vector d represents the change in distance and direction which an 
infinitesimal element of the medium experiences as a consequence of the applied stress. 

• In Cartesian coordinates with tensor notation, the last equation is given as: 

Ojd.j -f- djdi 

Hi = o 



5.5 Stress Tensor 

• The stress tensor, which is also called Cauchy stress tensor, is a rank-2 symmetric 29 
tensor used for transforming a normal vector to a surface to a traction vector acting on 
that surface, that is: 

T = <rn (250) 

where T is the traction vector, cr is the stress tensor and n is the normal vector. This is 
usually represented in tensor notation using Cartesian coordinates as: 


Tj (Tij Tlj 


(251) 


• The diagonal components of the stress tensor represent normal stresses while the off- 
diagonal components represent shear stresses. 


28 Some authors do not include the factor | in the definition of 7. 

29 In fact it is symmetric in many applications (e.g. in the flow of Newtonian fluids) but not all, as it can 
be asymmetric in some cases. We also choose to define it within the context of Cauchy stress law which is 
more relevant to the continuum mechanics; however it can be defined differently in other disciplines and 
in a more general form. 
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• Because the stress tensor is symmetric, in an nD space it possesses independent 

components instead of n 2 . Hence in a 3D space (which is the ordinary space for this 
tensor) it has six independent components. 

• In fluid dynamics, the stress tensor (or total stress tensor) is decomposed into two main 
parts: a viscous contribution part and a pressure contribution part. The viscous part may 
then be split into a normal stress and a shear stress while the pressure part may be split 
into a hydrostatic pressure and an extra pressure. 


5.6 Displacement Gradient Tensors 

• These are rank- 2 tensors which are denoted by E and A. They are defined in Cartesian 
coordinates using tensor notation as: 


E, 


v 


dxj 

dx 'j 


& 


Ay 


dp 

dxj 


(252) 


where x and x represent the Cartesian coordinates of an observed continuum particle at 
the present and past times respectively. These tensors may also be called deformation 
gradient tensors. 

• E quantifies the displacement of the particle at the present time relative to its position at 
the past time, while A quantifies its displacement at the past time relative to its position 
at the present time. 

• From their definitions, it is obvious that E and A are inverses of each other and hence: 



5.7 Finger Strain Tensor 


5.7 Finger Strain Tensor 

• This, which may also be called the left Cauchy-Green deformation tensor, is a rank-2 
tensor used in the fluid and continuum mechanics to describe the strain in a continuum 
object, e.g. fluid, in a series of time frames. It is defined as: 

B = E E t (254) 

which in Cartesian coordinates with tensor notation becomes: 

_ dxi dxj 
ij ~ dx' k dx' k 

where E is the first displacement gradient tensor as defined in § 5.6, the superscript T rep- 
resents matrix transposition, and the indexed x and x represent the Cartesian coordinates 
of an element of the continuum at the present and past times respectively. 


(255) 


5.8 Cauchy Strain Tensor 

• This, which may also be called the right Cauchy-Green deformation tensor, is the inverse 
of the Finger strain tensor and hence it is denoted by B T Consequently, it is defined as: 


B 1 = A r • A 

which in Cartesian coordinates with tensor notation becomes: 


(256) 


B -i = dx 'k 9x' k 
dxi dxj 


(257) 


where A is the second displacement gradient tensor as defined in § 5.6. 

• The Finger and Cauchy strain tensors may be labeled as “finite strain tensors” as 
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opposite to infinitesimal strain tensors. They are symmetric positive definite tensors; 
moreover they become the unity tensor when the change in the state of the object from 
the past to the present times consists of rotation and translation with no deformation. 

5.9 Velocity Gradient Tensor 

• This is a rank-2 tensor which is often used in fluid dynamics and rheology. As its name 
suggests, it is the gradient of the velocity vector v and hence it is given in Cartesian 
coordinates by: 

[Vv]jj = diVj (258) 

• The velocity gradient tensor in other coordinate systems can be obtained from the 
expressions of the gradient of vectors in these systems, as given, for instance, in § 4.3 and 
§ 4.4 for cylindrical and spherical coordinates. 

• The term “velocity gradient tensor” my also be used for the transpose of this tensor, i.e. 

(Vv) T . 

• The velocity gradient tensor is usually decomposed into a symmetric part which is the 
rate of strain tensor S (see § 5.10), and an anti-symmetric part which is the vorticity 
tensor S (see § 5.11), that is: 

Vv = S + S (259) 


5.10 Rate of Strain Tensor 


• This tensor, which is also called the rate of deformation tensor, is the symmetric part of 
the velocity gradient tensor and hence is given by: 30 


S = 


Vv + (Vv) 


(260) 


30 Some authors do not include the factor ! 2 in the definition of S and S and hence this factor is moved 
to the definition of Vv. Also these tensors are commonly denoted by 7 and u> respectively. 
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which, in tensor notation with Cartesian coordinates, is given by: 

c _ diVj + djVi 
'- , b — o 


( 261 ) 


• The rate of strain tensor is a quantitative measure of the local rate at which neighboring 
material elements of a deforming continuum move with respect to each other. 

• As a rank-2 symmetric tensor, it has n ^" 2 +1) independent components which is six in a 
3D space. 

• The rate of strain tensor is related to the infinitesimal strain tensor (refer to § 5.4) by: 


S = 


^7 

~dt 


(262) 


where t is time. Hence, the rate of strain tensor is normally denoted by 7 where the dot 
represents the temporal rate of change. 


5.11 Vorticity Tensor 

• This is the anti-symmetric part of the velocity gradient tensor and hence is given by: 


V v - (Vv) T 


(263) 


which, in tensor notation with Cartesian coordinates, is given by 

5 _ - djVi 

ij ~ 2 

• The vorticity tensor quantifies the local rate of rotation of a deforming continuum 
medium. 

• As a rank-2 anti-symmetric tensor, it has independent components which is three 


(264) 
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in a 3D space. These three components added to the six components of the rate of 
strain tensor give nine independent components which is the total number of independent 
components of their parent tensor Vv. 
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